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FOREWORD 


This  Technical  Report  was  prepared  in  the  Aeronautical 
Engineering  Department  of  the  University  of  Illinois,  Urbana, 
Illinois,  for  the  Fluid  Dynamics  Research  Branch,  Aeronautical 
Research  Laboratory,  Wright  Air  Development  Center,  on  Contract 
AF  33(6l6)-52.  Project  No.  3066,  "Gas  Turbine  Technology,"  Task 
No.  70154,  "Viscosity  Effects  in  Subsonic  Flow  Compressors  and 
Turbines."  with  Mr.  S.  Hasinger  acting  as  task  scientist.  This 
is  the  sixth  and  final  report  on  this  project. 

The  work  was  carried  out  under  the  direction  of  Dr.  T.  P. 
Torda,  Professor,  Project  Director,  who  is  responsible  for  the 
writing  of  this  report.  Several  workers  cooperated  on  this 
project  and  it  would  be  difficult  to  evaluate  the  contribution 
of  each  individual.  Instead.,  the  author  wishes  to  express  his 
sincere  thanks  to  his  co-workers  for  their  efforts  on  this 
project.  Appreciation  is  due  to  Dr.  H'.  H.  Hilton,  Assistant 
Professor,  Messrs.  P.  W.  Born,  F.  C.  Hall,  H.  A.  Hassan,  H.  R. 
Russell,  Research  Assistants,  and  to  the  undergraduate  computers 
who  worked  on  the  project. 
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ABSTRACT 


A  historical  review  of  the  various  turbomachine  theories  is  given*  The  Lorenz 
channel  theory  is  developed  for  the  analysis  of  viscous  flows  through  turboma chines. 
Both  the  infinitesimally  and  the  finitely  spaced  blades  are  considered.  First*  the 
incompressible*  viscous  and  then  the  compressible,  viscous  flows  are  treated.  The 
method  of  analysis  consists  of  an  iteration  procedure  in  which  the  first  step  is  the 
evaluation  of  an  equivalent  axially  symmetric  turbomachine  stage  with  infinites!* 
mally  spaced  blades*  Thus*  an  axially  symmetric  flow  field  is  used  and  the  stream¬ 
lines  are  determined  for  a  desired  pressure  rise  and  an  appropriately  chosen 
velocity  field  in  the  stags.  These  streamlines  form  the  infinitely  thin  streamline 
surfaoes*  l.e.  the  blades  of  the  axially  symmetric  case.  The  results  of  this  first 
step  in  the  analysis  are  used  as  the  first  terms  of  infinite  series  for  the  veloci¬ 
ties  and  pressure,  when  the  transition  to  finitely  spaced  blades  is  made.  These 
infinite  series  are  written  in  powers  of  the  tangential  coordinate  and  the  coeffi¬ 
cients  are  functions  of  the  radial  as  well  as  the  axial  coordinates.  The  fictitious 
force  field  of  the  analysis  of  infinitesimally  spaced  blades  is  eliminated  and  the 
complete  Navier-Stotes  equations  are  solved.  The  iteration  starts  from  a  number  of 
frozen  streamline  surfaces  obtained  from  the  analysis  of  the  axially  symmetric  case. 

There  are  a  number  of  boundary  conditions  which  must  be  satisfied  by  the 
solution  of  the  finitely  spaced  blade  analysis.  These  are>  1.  Along  the  frozen 
streamline  surfaces*  the  values  of  the  flow  variables  have  been  prescribed  to  be 
those  of  the  infinitesimally  spaced  blade  analysis.  2.  Along  the  new  and  corrected 
streamline  surfaces  forming  the  actual  blade  surfaces  in  the  finitely  spaced  blade 
analysis,  the  resultant  velocity  should  be  zero  (non-slip  condition).  This  latter 
boundary  condition  would  mate  the  already  non-linear  problem  more  non-linear  and 
render  the  solution  practically  Impossible.  Therefore*  these  conditions  have  been 
relaxed.  Since  in  modern  turbomaehlnes  the  blades  are  very  thin,  it  was  thought 
sufficiently  accurate  to  prescribe  that  the  non-slip  conditions  be  enforced  at  the 
chord  surfaces  of  the  new  blades  and  only  approximately  fulfilled  at  the  actual 
blade  surfaoes.  3.  The  conditions  of  closure  of  the  blade  surfaces  at  the  leading 
and  trailing  edges  of  the  blades  must  be  satisfied.  Numerical  illustrative  examples 
are  included.  Thirty-three  references  are  given. 
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INTRODUCTION 


Application  of  turbojet  engines  for  aircraft  propulsion 
intensified  the  efforts  to  develop  a  useful  method  for  the 
evaluation  and  design  of  turbomachines.  In  the  experimental 
field,  systematic  investigations  of  elements  and  also  of 
complete  engines  have  been  carried  out  in  order  to  determine 
the  behavior  of  airfoil  blades  in  cascade  arrangement.  The 
analytical  methods  employed  at  first  were  those  used  for 
steam  and  gas  turbines  and  were  of  a  semiempirical  nature 
using  average  velocity  distributions  between  blades  in  cas¬ 
cade  arrangement.  This  approach  may  be  called  the  momentum 
theory  and  it  gives  reasonably  good  results  for  the  first 
evaluation  of  turbomachines.  When  refinements  in  analysis 
were  required,  two  methods  emerged. 

The  first  method  of  approach,  the  direct  method, 
attempts  the  evaluation  of  the  flow  between  selected  air¬ 
foils  of  known  characteristics.  Conformal  mapping,  source 
and  sink  or  vortex  distribution  is  used  for  the  analysis  of 
the  flow  fields  and  pressure  distribution  around  the  blades. 
The  influence  of  the  cascade  arrangement  is  approximated  by 
empirical  factors  accounting  for  the  interaction  of  the  flow 
and  pressure  fields  around  adjacent  blades,  since  the  evalu¬ 
ation  of  this  effect  by  mathematical  methods  is  practically 
impossible.  Thus,  at  best,  only  an  approximation  of  the 

NOTE  -  This  technical  report  released  for  publication  in  November  1955. 
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interaction  between  blades  is  possible  with  this  method. 
However,  if  the  number  of  blades  is  small,  i.e.  the  inter¬ 
action  between  blades  is  small,  the  direct  method  can  yield 
useful  data  for  turbomachine  analysis. 

The  second  method  of  approach  is  the  channel  theory,  or 
indirect  method.  The  passages  between  blades  are  treated  as 
channels  and  the  shapes  of  their  boundaries,  i.e.  the  shape 
of  the  blades,  are  determined  from  the  desired  distribution 
of  the  flow  variables  (the  velocity  and  pressure  distribu¬ 
tions)  through  the  individual  stages.  It  is  of  interest  to 
note  that  the  channel  theory  was  developed  first.  When  the 
circulation  theory  was  introduced  for  analysis  of  cascades, 
the  channel  theory  was  forgotten,  probably  due  to  the  fact 
that  in  its  early  form  of  development  it  was  suitable  only 
for  the  treatment  of  a  turbomachine  stage  with  infinite 
number  of  blades  and  only  approximate  corrections  could  be 
made  to  account  for  a  finite  number  of  blades. 

In  the  following  brief  historical  review  is  given  the 
development  of  analytical  methods  for  the  evaluation  and 
design  of  turbomachines .  No  attempt  is  made  to  present  a 
complete  review,  but  rather  the  early  development  is  treated 
in  more  detail.  This  is  done  because  the  more  recent  work 
is  well  known  by  men  working  in  this  field.  Some  represent¬ 
ative  references  are  included  which  seem  to  have  been  of 
importance  in  the  development  of  the  channel  theory.  This 
list  should  not  be  considered  as  complete  by  any  means. 
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The  development  of  the  theories  of  turbomachinery  can 

be  traced  to  analyses  of  water  turbines  and  ship  propellers 

which  were  the  first  practical  applications  of  rotating 

flow-machinery.  The  first  important  step  in  the  formulation 

of  a  theory  was  made  by  L.  Euler  in  about  1740  with  the 

derivation  of  the  equations  of  motion  for  a  nonviscous, 

homogeneous,  incompressible  fluid.  Somewhat  later  he 

discussed  the  theory  of  water  turbines  proposed  by  Segner  in 
1* 

1750  .  He  also  gave  a  somewhat  more  general  treatment  of 

the  subject  in  subsequent  papers2*^.  F.  Prasil  applied 
Euler *s  equations  of  motion  to  the  investigation  of  water 
turbines  in  cylindrical  coordinates  in  1903^.  Inspired  by 
Prasil's  work,  H.  Lorenz  developed  the  theory  of  water  tur¬ 
bines  and  ship  propellers  and  published  his  work  in  several 
papers  in  1905  and  1906.5>^»7»^*9»10*1]-  The  results  are 
given  in  his  famous  book,  A  New  Theory  and  Calculations  of 
Turbine  Wheels,  first  published  in  1906  •  Lorenz  based  his 

investigations  on  the  analysis  of  an  axially  symmetric  flow 
field  in  which  the  streamline  surfaces  correspond  to  an 
infinite  number  of  blades,  hence  the  axial  symmetry  of  the 
flow.  He  introduced  the  concept  of  "forced  accelerations" 
(Zwangsbeschleunigung) ,  i.e.  a  force  field  per  unit  of  mass 
and  represented  the  effect  of  the  infinite  number  of  blades 
on  the  working  fluid  by  these  accelerations.  Thus,  the 

^Superscripts  denote  References ' 
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"jump"  in  pressure  distribution  on  either  side  of  the  infin¬ 
itesimally  spaced  "blades"  was  represented  by  the  force 
field*  W.  Bauer sf eld  showed  that  the  force  field  must  be 
perpendicular  to  the  streamlines  for  a  nonviscous  fluid^. 
The  concept  of  the  Lorenz  "forced  accelerations"  has  been 
criticized  severely  by  Prasil  and  von  Mises,  but  Stodola^, 
Bauersfeld^,  and  Loewy^  have  successfully  defended  Lorenz 
and  have  shown  that  the  attacks  on  his  theory  were  based  on 
misunderstandings*  In  1909*  R.  von  Mises  published  his 
"Theory  of  Water  Wheels"1^  and  derived  equations  similar  to 
those  of  Lorenz.  In  1912,  Bauersfeld  applied  the  Lorenz 
theory  to  the  analysis  and  design  of  a  Francis  wheel  and, 
thus,  established  its  usefulness  for  practical  applications. 

In  the  meantime,  Kutta  applied  the  Joukowsky  concept  of 
circulation  to  the  analysis  of  a  straight  lattice  of  airfoil 

blades1^  and  in  1913,  H.  Foettinger  introduced  the  use  of 

20 

circulation  in  turbine  theory  .  Many  workers,  among  them 
Bauersfeld,  Koenig,  and  Spannhake,  applied  the  circulation 
theory  to  plane  stationary  lattices  of  airfoils.  However, 
the  momentum  theory  which  is  given  in  Stodola's  "Steam  and 
Gas  Turbines"  has  become  the  standard  method  for  the 
evaluation  of  turbomachinery*  Since  this  theory  is  based 
on  average  velocities  over  cross  sections,  a  first  approxi¬ 
mation  of  design  and  output  of  the  machines  is  possible 
only  with  this  method* 
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As  far  as  the  author  is  aware  of  developments,  H.  J. 
Reissner  in  the  late  1940 »s  had  revived  the  Lorenz  theory  to 
the  analysis  of  turbomachinery  and  stationary  lattices.  Al¬ 
though  Lorenz  had  already  introduced  corrections  for  a  finite 
number  of  blades  in  his  analysis,  Reissner  modified  the  method 
and  introduced  the  transition  from  the  infinitesimally  spaced 
blades  to  those  of  finite  spacing  by  "freezing"  a  certain  num¬ 
ber  of  streamline  surfaces  obtained  from  the  axially  symmetric 
analysis  and  calculating  the  flow  between  them  by  a  perturba¬ 
tion  method.  Reissner  and  his  collaborators  published  several 
papers22*2^ *2^  on  the  flows  between  finitely  spaced  lattices 
for  incompressible  and  for  compressible  non- viscous  flows. 

In  1949,  the  author  suggested  a  method  of  application  of  the 
Lorenz-Reissner  theory  to  the  analysis  of  turbomachines  with 
viscous  flows2^.  This  approach  has  been  applied  by  Reissner 
and  his  co-workers  to  viscous  flows  between  stationary  lat- 
tices  .  However,  all  of  Reissner’s  work  has  been  done  with 
a  one-step  perturbation  only  and  certain  boundary  conditions 
(the  closure  of  the  blade  surfaces  at  the  leading  edge  and 
the  trailing  edge)  could  not  be  satisfied. 

Work  at  the  NACA  by  C,  Stanitz,  C.  H.  tfu  and  others 
should  also  be  mentioned,  as  well  as  analyses  by  F.  E.  Marble, 
V.  D.  Rannie,  and  others  at  California  Institute  of  Technology. 
In  nearly  all  work  done  to  date  the  direct  method  of  approach 
has  been  used,  or  a  mixture  of  direct  and  indirect  methods, 
and  in  most  cases  axial  symmetry  of  the  flow  was  assumed. 

Some  numerical  methods  for  correction  for  finite  blade  spacing 
exist,  but  a  general  treatment  for  finite  blade  spacing  has 
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not  boon  given.  As  far  as  the  author  knows ,  the  first  com¬ 
plete  analysis  of  viscous  flows  through  turbomachines  was 
done  on  this  contract  for  the  case  of  infinitesimally  spaced 
blades^  and  also  for  the  case  of  finitely  spaced  blades^. 

The  analysis  developed  under  this  contract  for  the  eval¬ 
uation  and  design  of  turbomachines  is  based  on  the  Lorenz- 
Reissner  channel  theory.  The  method  consists  of  an  iteration 
procedure  in  which  the  first  step  is  the  evaluation  of  an 
equivalent  axially  symmetric  turbomachine  stage  with  infini¬ 
tesimally  spaced  blades.  Thus,  an-  axially  symmetric  flow 
field  is  used  and  the  streamlines  are  determined  for  a  desired 
pressure  rise  and  an  appropriately  chosen  velocity  field  in 
the  stage.  These  streamlines  form  the  infinitely  thin  stream¬ 
line  surfaces,  i.e.  the  blades  of  the  axially  symmetric  case. 
The  results  of  this  first  step  in  the  analysis  are  used  as 
the  first  terms  in  infinite  series  for  the  velocities  and 
pressure.  These  infinite  series  are  written  in  powers  of  the 
tangential  coordinate  and  the  coefficients  are  functions  of 
the  radial  as  well  as  the  axial  coordinates. 

The  basis  for  the  iteration  is  the  fact  that  for  the 
axially  symmetric  case  a  force  field  per  unit  of  volume  is 
assumed  which  represents  the  "jumps"  of  the  pressure  on  either 
side  of  the  streamline  surfaces,  i.e.  the  force  field  repre¬ 
sents  the  action  of  the  "blades"  on  the  working  fluid.  This 
force  field  is  eliminated,  in  the  analysis  of  the  flow  through 
the  stage  with  a  finite  number  of  blades,  and  the  additional 
terms  of  the  series  correct  the  streamlines  and  evaluate  the 
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deviations  from  the  axially  symmetric  streamlines.  The 

corrected  streamline  surfaces  form  the  upper  and  lower  faces 
of  the  actual  finitely  spaced  blades.  Along  with  the  correc¬ 
tion  of  the  streamlines,  the  velocity  and  pressure  distribu¬ 
tions  of  the  axially  symmetric  case  are  corrected  as  well  and 
the  flow  field  in  the  stage  with  finitely  spaced  blades  is 
evaluated.  The  method  of  transition  from  the  axially  symmet¬ 
ric  to  the  axially  non-symmetric  case  consists  of  "freezing" 
a  certain  number  of  streamline  surfaces  of  the  axially  symmet¬ 
ric  case  and  assuming  that  these  frozen  streamline  surfaces 
will  lie  somewhere  near  the  center  of  the  channel  between  two 
blades.  The  freezing  means  that  along  these  surfaces  the 
flow  in  the  finite  case  is  prescribed  to  be  the  same  as  in 
the  infinitesimal  case. 

In  the  analysis  of  the  axially  symmetric  case,  the 
Navier-Stokes  momentum  equations  (including' the  terms  of  the 
Lorenz  force  field  for  a  viscous  fluid)  and  the  continuity 
equation  are  used  for  the  determination  of  the  flow  variables 
for  prescribed  inflow  and  exit  conditions.  The  type  of  the 
resultant  velocity  profile  in  the  channel  between  the  infin¬ 
itesimally  spaced  blades  is  assumed  and  the  Navier-Stokes 
equations  integrated.  Then  the  streamline  equations  are 
integrated  along  the  stage.  The  pressure  and  velocity  dis¬ 
tributions  along  the  blade  channels  can  be  plotted  once  the 
intake  and  exit  conditions  are  satisfied  and  the  desired 
variation  of  the  flow  variables  achieved  by  judicious 
choice  of  the  arbitrary  functions  of  integration. 
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The  same  equations  are  used  for  the  analysis  of  finitely 
spaced  blades,  except  that  the  Lorenz  force  field  is  elimin¬ 
ated  from  the  Navier-Stokes  equations  and  infinite  series  are 
introduced  for  the  flow  variables  (the  three  velocity  compon¬ 
ents  and  the  pressure ).  By  substituting  the  series  into  the 
Navier-Stokes  and  the  continuity  equations,  recurrence  formu¬ 
lae  are  obtained  for  the  coefficients  of  the  series.  It  was 
mentioned  that  the  first  terms  of  the  series  (  u0  ,  irg  , 
and  -j&o  )  are  those  which  were  obtained  from  the  analysis  of 
the  infinitesimally  spaced  blades.  The  series  of  two  out  of 
three  of  the  velocity  variables  contain  arbitrary  functions 
as  the  second  terms  of  the  series.  Thus,  u.,  and  either  v? 
or  w,  have  to  be  determined  from  boundary  conditions.  Once 
the  coefficients  of  the  series  are  determined,  the  streamline 
equations  may  be  integrated.  They  yield  the  corrected  stream¬ 
line  surfaces  forming  the  upper  and  lower  faces  of  the  actual 
blades. 

There  are  a  number  of  boundary  conditions  which  must  be 
satisfied  by  the  solution  of  the  finitely  spaced  blade  anal¬ 
ysis.  These  are:  1.  Along  the  frozen  streamline  surfaces, 
the  values  of  the  flow  variables  have  been  prescribed  to  be 
those  of  the  infinitesimally  spaced  blade  analysis.  2.  Along 
the  new  and  corrected  streamline  surfaces  forming  the  actual 
blade  surfaces  in  the  finitely  spaced  blade  analysis,  the 
resultant  velocity  should  be  zero  (non-slip  condition).  This 
latter  boundary  condition  would  make  the  already  non-linear 
problem  more  non-linear  and  render  the  solution  practically 
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impossible*  Therefore,  these  conditions  have  been  relaxed* 
Since  in  modern  turbomachines  the  blades  are  very  thin,  it 
was  thought  sufficiently  accurate  to  prescribe  that  the  non¬ 
slip  conditions  be  enforced  at  the  chord  surfaces  of  the  new 
blades  and  only  approximately  fulfilled  at  the  actual  blade 
surfaces*  This  means  that  the  resultant  velocity  at  the  blade 
surfaces  will  not  be  zero  but  have  a  very  small  value.  This 
relaxation  of  the  non-slip  conditions  makes  the  boundary  con¬ 
ditions  linear*  3*  The  conditions  of  closure  of  the  blade 
surfaces  at  the  leading  and  trailing  edges  of  the  blades  must 
be  satisfied*  The  closure  at  the  leading  edges  can  be  enforced 
by  judicious  choice  of  arbitrary  integration  constants  and  the 
closure  at  the  trailing  edges  by  the  correct  choice  of  the 
arbitrary  functions  LL(,  and  (or  VJ',  ).  At  first,  the 

analysis  of  flows  through  turbomachine  stages  with  infinites¬ 
imally  spaced  blades  will  be  presented*  The  basic  idea  for 
such  an  analysis  is  given  in  Ref*  25  which  made  the  Inclusion 
of  viscous  effects  possible*  This  method  was  applied  first 
to  laminar.  Incompressible  flows2"'7  and  somewhat  generalized 
later  along  the  lines  suggested  by  author  to  include  laminar 
as  well  as  turbulent  compressible  or  incompressible  flows2 
Numerical  Illustrative  examples  are  also  included* 

Then,  the  extension  of  the  analysis  to  turbomachine 
stages  with  finite  number  of  blades  is  given  for  incom¬ 
pressible2®*^®  and  compressible^  flows  and  the  difficulties 
of  considering  turbulence  are  discussed* 
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ANALYSIS 


I.  VISCOUS  FLOWS  THROUGH  TURBOMACHINE 
STAGES  WITH  LARGE  NUMBER  OF  BLADES 
INFINITESIMAL  BLADE  SPACING. 

1.1.  METHOD  OF  SOLUTION. 

In  considering  viscous  flows  through  turbomachine 
stages  with  infinitesimally  spaced  blades,  the  Lorenz  force 
field  is  generalised  to  include  the  action  of  shear  forces. 
The  force  field,  as  conceived  by  Lorenz,  accounted  for  the 
action  of  the  "blades"  on  the  fluid.  This  means  that  only 
the  difference  of  pressure  on  either  side  of  the  streamline 
surfaces  was  represented  by  Inclusion  of  force  field  terms 
on  the  momentum  equations. 

The  assumption  of  infinitesimally  spaced  blades  is 
equivalent  to  axial  symmetry  of  the  flow  and  it  also  allows 
the  use  of  the  average  resultant  velocity  between  these 
blade  surfaces  instead  of  the  actual  profile  of  the  result- 
ant  velocity. 

The  analysis,  as  presented,  is  valid  for  rotating  coor¬ 
dinate  systems  (rotating  stage).  For  a  stator  stage,  the 
Coriolis  force  terms  are  equal  to  zero. 

The  Navier-Stokes  equations  for  axially  symmetric  flows 
(  y|r  *  0)  written  in  cylindrical  coordinates  and  including 


WADC  TR  55-231 


10 


? 


I  f 


the  Coriolis  forces  and  the  force  field  terms  are 


u!^  +  nr&t  ~-vlir  -  2irai-cj*r  -  -7*  fr  **f 


UL 


^  +iufc)  =  -f> 


(1.1.1) 


U&+vrtoL  *-?d£  *-f= 


where 

i jL,ir,iLr 

*> 

«o 

i 


are  the  cylindrical  coordinates 

are  the  velocity  components  in  the  r,<¥ 
and  z  directions,  respectively.  These 
velocity  components  are  measured  relative 
to  the  rotating  coordinate  system 

is  the  pressure 

is  the  density 

is  the  constant  angular  velocity 

is  the  force  field  component  per  unit  of 
mass 


subscripts  r,  9  and  z  denote  directions 
Fig.  1  gives  the  notation  used. 

The  force  field  accounts  for  the  effects  of  the  dilata¬ 
tion  and  the  shear  stresses.  Its  components,  ff  ,  f  <p  and 
f  z  ,  are  inclined  to  the  streamline  surfaces  and  each  may 
be  written  as  the  vector  sum  of  two  components,  ■** 
etc.,  such  that  and  are  force  components  per 
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unit  of  mass  perpendicular  to  the  streamline  surfaces ,  and 
i-  r ,  and  'ft  z  are  force  components  per  unit  of  mass 

tangent  to  the  streamline  surfaces. 

The  magnitude  o£  the  resultant  velocity,  relative  to 
the  rotating  coordinate  system,  is 

(1.1.2) 

As  was  mentioned |  instead  of  the  actual  value  of  the 
resultant  velocity  V  ,  its  average  will  be  used  in  the 
analysis.  This  average  may  be  determined  depending  on  the 
laminar  or  turbulent  nature  of  the  flow.  The  procedure, 
though  well  known,  will  be  explained  later  in  the  chapter 
on  "Numerical  Examples".  Therefore,  from  here  on,  V  will 
denote  the  average  of  the  resultant  relative  velocity  at 
any  point  between  the  streamline  surfaces.  The  resultant 
of  the  force  field  due  to  shear  forces  'A.  and  its  com¬ 
ponents  satisfy  the  relationship 

■fu  \/  -  Vft_r  XL  +  V"  *  XLf  (1.1.3  ) 

and  the  y  -  components  of  ft  satisfy 

ar-u.  ^ir  f  (^ur  =  O  (1.1.4) 

since  they  are  perpendicular  to  V  • 

Multiplication  of  the  Navier-Stokes  equations  (1.1.1) 
by  u.  ,  ‘Xr  and  -ur  ,  respectively,  and  summation  of  the  three 
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component  equations ,  yields 


T)K  _  n 

Hx  =  %  V 


where  t>t 


P  _  Z 
D-t  'LL  Sr 


U r 


A 


and  \  s  -%r  “  -&4aJt"  +  P 


(1.1.5) 


(1.1.6) 


with  P  ,  the  enthalpy  function, 

p.yv 

Eq.  (1.1.5)  may  be  re-written  as 

u4t“  +  7xr££  =  "fu  V 


(1.1.7) 


(1.1.6) 


and  eq.  (1.1.3)  is  a  Lagrange  equation  which  may  be  solved 
by  finding  a  solution  of  the  following  system  of  ordinary 
differential  equations,  the  subsidiary  equations  of  (1.1.3) 


-  sLz  -  d.  K. 

ix.  us  vCv 


(1.1.9) 


The  complete  Integral  of  eq.  (1.1.3)  is 

(,r,x.^=0  (1.1.10) 

where 

t  (x\1,  N)-oL  and  •?  (r.z,  O  =  A  (1.1.11) 

are  two  independent  solutions  of  eqs.  (1.1.9),  the  charac¬ 
teristics  of  the  system,  oL  and  A  are  constants  of 


WADC  TR  55-231  13 


j 


V 


integration,  and  V  is  an  arbitrary  function. 

The  first  equation  in  (1.1.9) ,  i.e. 

dir  _  d.  z. 

~Ul  "ur 


(1.1.12) 


ia  the  equation  of  the  r-z.  trace  of  the  streamlines.  Let 
the  integral  of  eq.  (1.1.12)  be 

5  (t.z.'J*  <t  (1.1.13) 


where  <r  ia  a  constant  of  Integration  and  to  each  value  of 
(jr  corresponds  an  r-x  trace  of  a  streamline. 

Eq.  (1.1.15)  (and,  therefore,  eq.  (1.1.8))  must  be 
integrated  along  streamlines,  since  they  represent  the 
resultant  of  the  component  momentum  equations.  Let  in 
(1.1.11) 


(1.1.14) 


and 


•*-/ 


Ijjr  CL  2. 


Let 

R(r,xY»  j*' d.'r  and  T.ir.x)^  j 

then,  from  (1.1.14)  and  (1.1.15) 


(1.1.15) 


^dj.  (1.1.16) 


A”  ft  t cL 


(1.1.17) 
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and 


X  *  Z  (r.z')  +■  P  (1.1,18) 

However,  through  (1.1.10),  (1.1.11),  (1.1.12)  and 
(1.1.13),  a  functional  relationship  exists  between  ec  , 
and  $  •  Let 

o csG,(<r)  and  ^=0^(0-^  (1.1.19) 

then,  eqs.  (1.1.17)  and  (1.1.18)  may  be  written  as 

X  =  R  Cr.x')  **-0,(0  (1.1.20) 

and 

X  *  Z  ( r.iO  +  ©i,  (^"i)  (1.1.21) 

since  the  integration  must  be  carried  out  along  streamlines 
(here,  along  r-z.  traces  of  the  streamlines). 

For  the  solution  of  (1.1.18),  the  values  of  X  must  be 
given  along  an  initial  curve,  C  (re (ffOjZoGr))  ,  e.g.  the 
leading  edge  of  the  compressor  blade.  This  means  that 

X  =  X0  (Ve)z0')  (1.1.22) 

must  be  known  and  the  only  condition  on  C(YClx0}  is  that 
it  should  not  be  a  streamline  trace  (characteristic  line). 

The  initial  value  of  X  defines  the  functions  (o') 
and  (O’)  •  If  f  is  eliminated  between  eqs.  (1.1.20) 
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and  (1.1.21),  the  solution  X*X('»\x>)  is  obtained  which 
satisfies  the  differential  equation  (1.1.8)  and  the  initial 
conditions  (1.1.22). 

The  solution  of  eq.  (1.1.8),  as  described,  is  straight¬ 
forward,  if  ix  ,  *ur  and  '■ft  V  are  known  functions  of  v  , 

Z  and  X  ,  and  if  -u.  ,  us  and  'ft.  V  and  their  first 
derivatives  are  continuous  in  the  T  -  Z  -  X  space  of 
interest.  As  a  matter  of  fact,  the  channel  theory  requires 
the  determination  of  the  blade  shapes  from  desired  velocity 
and  pressure  distribution  through  the  stage. 

In  what  follows,  a  somewhat  more  general  treatment  of 
the  problem  is  given.  Instead  of  prescribing  "U.  , 

hr(r,7L)  and  vp(r,z')  (and,  therefore  K(r,  Xs)  ), 
the  types  of  R  (>,*.')  and  Z  (r,x)  »  see  eqs.  (1.1.16), 

will  be  prescribed  and  these  functions,  together  with  the 
initial  conditions,  will  determine  the  flow  parameters  as 
well  as  the  blade  shapes  throughout  the  stage.  Assumption 
of  separable  variables  in  R  and  Z.  greatly  simplifies  the 
solution  and  yields  good  flow  parameter  distribution  and 
blade  shapes. 

For  the  complete  solution  of  the  problem,  the  continu¬ 
ity  equation 


must  be  satisfied  also,  and  either  an  equation  of  state. 
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•p  —  '-p(e)T')  ,  or  a  -p-P  relationship  must  be  given 

in  form  of  a  change  of  state  equation*  '-p  =  . 

1.2.  SPECIAL  SOLUTIONS. 

In  the  following,  some  special  solutions  are  given. 

It  will  be  seen  that  reasonably  general  cases  are  obtained 
even  though  simplifying  assumptions  are  made. 

1.2.  a.  ANALYSIS  OF  LAMINAR  INCOMPRESSIBLE  FLOWS  WITH 

CONSTANT  AXIAL  VELOCITY  COMPONENT. 

It  is  usual  to  design  axial  flow  compressors  with  a 
constant  axial  velocity,  i.e. 

'ur  -  zltc  (1.2. a. 1) 

For  incompressible  flow,  the  density  6  ■  const.  Since  the 
flow  is  laminar,  one  may  assume  that  the  velocity  distribu¬ 
tion  in  each  channel  between  the  blades  is  a  non-symmetric 
parabola  (Fig.  2)  of  the  non-dimensional  form* 

y)  =•  (/-  $z)/(l  -,2a.§  +CLa)  and  |al<l  (I.2.a.2) 

where 


y)  =  V/SAno*  %  *  y)/b 

a  =  e/b  b  =  d/a. 


(I. 2. a. 3) 


and  e,  is  the  eccentricity. 


♦the  concept  of  parabolic velocity  distribution  in  the 
channel  was  first  introduced  in  Ref.  6. 
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It  should  be  noted  that  the  parabola  vised  here  (non- 
symmetric  with  respect  to  the  resultant  velocity,  V  )  is 
representative  only.  Other,  higher  order  polynomials  may  be 
used  instead  if  turbulent  flows  are  investigated  and  a  flat 
velocity  distribution  exists  between  the  blades.  However, 
the  choice  of  the  type  of  velocity  distribution  is  signifi¬ 
cant  only  for  the  results  of  the  case  of  infinitesimal  blade 
spacing,  since  this  velocity  distribution  is  corrected  by 
the  series  when  the  finite  blade  spacing  is  treated.  There¬ 
fore,  only  the  convergence  of  the  series  will  be  affected  by 
the  original  choice  of  the  curve  which  represents  the  result¬ 
ant  velocity  profile  between  blades. 

For  laminar  flow  the  shear  stress  equals 

T  -  y-v  {[dv/a-n3x=  4  +C^v/d(— a\]Tw  =  -  (1.2. a.4) 

where  -n.  is  the  direction  normal  to  the  streamlines 
d  is  the  blade  spacing 
/<  is  the  viscosity  (taken  as  constant). 

The  blade  spacing  (Fig.  1)  is 

oL=  r  (I.2.a.5) 

and  is  the  pitch  angle. 

The  resultant  shear  force  per  unit  volume  may  now  be 
calculated  from  eqs.  (I. 2. a. 4)  and  (I.2.a.2) 


Id 
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(I. 2. a. 6) 


where 


Vrnax  ”  3  \4.ve  ML 


and 


and 


(I.2.a*7) 


(I.2*a.$) 


(1,2. a. 9) 


for  a.  *  0  :  ML  =  I 

G,=  I  *.  ML  -  41 


(1.2. a. 10) 


The  analysis  which  follows  is  valid  for  any  prescribed' 
velocity  profile,  since  only  the  value  of  the  constant  ML 
varies  for  various  velocity  profiles* 

Eq.  (1.1*3)  with  (l*2.a.6)  becomes 

urf2  -GL^yP" ^ -  (I*2.a*ll)‘ 

For  incompressible  flow,  and  with  the  assumption  of  con¬ 
stant  axial  velocity  through  the  blade,  eq*  (I.2*a*l),  the 
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continuity  equation  (1*1.23)  reduces  to 
&r  '  O 


(I. 2. a. 12) 


with  the  solution 


U*Z,(x)/r 


(I. 2. a. 13) 


where  Z,  (z.)  is  an  arbitrary  function  of  z.  to  be 
determined  by  the  boundary  conditions  of  the  problem* 
Gq*  (1*1*5)  will  be  for  this  case 


D-t  "  pslb Vz-oA)4 


(I. 2. a. 14) 


If  the  quantity  CL  ,  eq*.  (I.2.a*3),  is  taken  as  con¬ 
stant,  then  all  the  points  of  maximum  velocity  Vty >o.x  in 
the  channel  lie  along  a  radius  at  any  particular  z  •  Sub¬ 
stitution  far  bx  from  eq.  (I. 2. a. 3)  and  use  of  eqs. 

(I. 2. a. 5)  and  (I*2*a*3)  reduces  eq*  (I*2.a*14)  to 

ti 


*  UTo'&V  ” 


(I. 2. a. 15) 


with 


(I.2*a*l6) 


Eq.  (I*2*a*15)  is  the  Lagrange  partial  differential 
equation  for  this  case  and  its  subsidiary  equations  are 

4e  «  4k  - 


(I. 2. a. 17) 
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One  Integral  of  (I. 2. a. 17)  la 


ol  z  +•  C , 


inhere  C,  la  an  arbitrary  constant. 
Let 


ra  "  dz  and  ©  =  © (<\x> 


(I. 2. a. 18) 


(I. 2. a. 19) 


then  eq.  (I. 2. a. 18)  becomes 

A  =4^4'Ci  (I.2.a.20) 

Eliminating  A  and  with  help  of  eqs. 

(I. 2. a. 19)  and  (.2. a. 20),  the  differential  equation 
(I. 2. a. 15)  simplifies  to 

\L  —  0  (I. 2. a. 21) 


Since  u.  &  o  ,  the  solution  is 
e  =  Zj.(7.) 


(I. 2. a. 22) 


and  2  *  (z)  is  an  arbitrary  function  of  z  along  each 
r  -  z  streamline  trace.  In  general,  ~  Z^.  ( . 
With  the  above  results,  the  solution  of  the  problem  is 
completed.  Thus 


(I.2.a.23) 
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(I. 2. a. 24) 


P  *=  +  vre  +  oj1  rx/2.  +■  C.,*J 


XT 


1* 


Tz 


(I. 2. a. 25) 


The  velocity  and  pressure  fields  in  the  compressor  are 
completely  described  by  the  foregoing  equations.  The  traces 
of  the  streamlines  may  be  determined  from 


sL£  -  rd<p  _  sLL. 
U  "  ir  Mr 


(I. 2. a. 26) 


The  integration  of  pairs  of  differential  equations  will  then 
yield  the  streamline  traces  in  any  plane. 

For  any  given  stage  of  the  compressor,  the  angular 
velocity  co  ,  the  magnitude  and  direction  of  the  entrance 
velocity,  the  entrance  pressure  and  hub  and  shroud  radii  at 
the  entrance  are  prescribed.  These  known  quantities  are 
used  to  determine  the  arbitrary  functions  Z, (z)  and 
0,  (r,x)  and  the  constant  C,  . 

It  should  be  remembered  that  since  these  latter  func¬ 
tions  are  arbitrary  it  is  permissible  to  assume  their  type. 
Such  assumptions  in  no  way  detract  from  the  uniqueness  of 
the  solution,  since  for  any  given  combination  of  2,9,  and 
C,  one  and  only  one  blade  profile  is  found.  Different 
velocity  and  pressure  fields  will  naturally  require 
different  blade  shapes. 
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C,  and  the  Initial  values  of  Z ,  (o)  and  Q  (o) 
are  determined  from  the  given  entrance  conditions.  The 
variation  of  Z,  and  0  from  entrance  to  exit  has  to  be 
chosen  to  conform  with  the  physical  limitations  on  pressure 
distribution  through  the  stage  and  the  curvature  of  the 
infinitesimal  "blades".  Too  rapid  a  rise  of  pressure 
corresponds  to  extreme  curvature  of  the  "blade"  surface  at 
the  leading  edge  which  would  increase  viscous  losses  and 
facilitate  flow-separation.  The  functions  "Z ,  and  0  also 
determine  the  total  twist,  (-0  ,  of  the  "blade"  surfaces. 

Excessive  values  of  PU)  also  increase  the  viscous  losses 
due  to  the  increase  in  "blade"  surface  area.  There  is  no 
method  available  for  the  a  priori  determination  of  the  types 
of  functions  to  be  used  for  Z  ,  and*©  .  Since  the  infin¬ 
itesimal  "blade"  surfaces  will  serve  as  chord-surfaces  of 
the  finitely  spaced  blades,  the  choice  of  Z ,  and  0  func¬ 
tions  may  be  made  to  obtain  known  efficient  chord  curvatures 
for  the  finitely  spaced  blades. 

I.2.b.  NUMERICAL  EXAMPLE  FOR  LAMINAR  INCOMPRESSIBLE  FLOWS 
WITH  CONSTANT  AXIAL  VELOCITY  COMPONENT 

The  aim  of  the  numerical  work  was  to  determine '  under 
which  conditions  a  maximum  pressure  rise  could  be  achieved 
through  one  rotating  stage  of  infinitesimally  spaced  blades 
(stream-surfaces).  At  the  same  time  large  twist  and  large 
curvatures  of  these  surfaces  were  to  be  avoided.  A  number 
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of  numerical  examples  were  worked;  however,  only  one  is 
presented  in  this  report  showing  trend  and  limitations  of 
the  design. 

There  are  two  arbitrary  functions  and  Q  (V>z)  , 

and  one  arbitrary  constant  C,  in  the  analysis.  C,  may 
be  determined  from  eq.  (I. 2. a. 24)  with  the  intake  conditions. 
Thus, 


c,=  $ 


AZjo) 

Zd~o 


(I.2.b.l) 


where  subscript  i.  denotes  values  at  the  intake  of  the 
rotating  stage • 

The  choice  of  the  arbitrary  functions,  r2^\(x)  and 
Q  (r^T.)  ,  determines  the  blade  shape. 

Eq.  (1.2. a.  13)  sheets  that  the  variation  of  the  radial 
component,  XL  ,  of  the  resultant  velocity,  V  ,  depends  on 
the  arbitrary  function  Z,  .  Through  the  equation  of  the 
streamlines,  eq.  (I.2.a.26),  2  ,  also  determines  the  r  -Z 
variation  of  shroud  and  hub.  The  requirement  of  small 
curvature  of  the  stream  surfaces  restricts  the  choice  of 
the  type  of  functions  used  for  Z ,  • 

The  pressure  rise  through  the  rotating  stage  along 
each  streamline  is  due  to  dynamic  recovery  and  centrifugal 
action. 

The  pressure  rise  due  to  dynamic  recovery  is  a  maximum 
if  the  resultant  velocity  at  the  exit,  Ve  *  *®r°»  (Sub- 
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script  e  denotes  exit  conditions.)  Since,  however,  the 
axial  component  of  the  relative  velocity,  10%  ,  is  taken 
as  constant  and  since  the  influence  of  the  radial  velocity 
component  is  small,  the  tangential  velocity  component  at 


the  exit, 

,  was  taken  as  zero.  Thus 

^  ♦  [ 

^rT"r 

(I.2.b.2) 

where  JL  is  the  value  of  z  at  exit.  With  the  condition  of 
ire  ,  the  stream-surface  at  the  exit  is  tangential  to  the 
T  -  2.  plane. 

The  pressure  rise  due  to  centrifugal  action  is  directly 
dependent  on  the  ratio  of  the  exit  and  intake  radii,  i.e., 
Te/vi  ,  and  on  the  angular  velocity  ^  •  Use  of 
eq.  (I.2.b.l)  in  (I.2.b.2)  yields 

P  _  fp;  ,  (*ir{'L  -  (co*- -  AZ,. 

<°  [  (°  <2.  J  X  (I.2.b.3) 

The  angular  velocity  is  prescribed  for  any  particular 
design  problem.  Then,  the  pressure  rise  due  to  centrifugal 
action  is  a  function  of  the  exit  radius,  re  •  If  this 
radius  is  chosen,  e.g.,  the  exit  hub-radius  might  be  pre¬ 
scribed  as  an  allowable  maximum,  then  the  maximum  pressure 
is  defined.  The  shape  of  each  streamline  in  the  r  -  z. 
plane  is,  thus,  completely  determined  by  prescribed  values 
of  the  respective  and  by  a  compromise  between  the  value 

of  a  maximum  pressure  rise  and  a  reasonable  hub  exit  radius. 
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and  by  the  requirement  of  gentle  variation  from  intake  to 
exit*  The  choice  of  Z ,  yielding  a  straight  conical  shroud 
and  Slightly  conical  hub  shape  is  thought  to  be  appropriate, 
and  this  is  the  reason  for  the  choice  of  the  function 


(I.2.b*4) 


where  subscript  £  denotes  shroud*  Note:  For  the  case  of 
straight  conical  shroud  S^sAz  ■  const*  *  Ms  • 

Once  the  r  -  Z  variation  of  the  streamlines  is  given, 
the  twist  of  the  stream-surface  is  defined  by  ©  ,  since 

=  .  In  order  to  reduce  the  length  of  the 

stream-surfaces,  and  thus,  reduce  the  magnitude  of  the  vis¬ 
cous  dissipation,  the  twist  angle  should  be  small*  This 
requirement  and  also  the  necessity  of  avoiding  sharp  curva¬ 
tures  of  the  stream-surfaces  allow  the  choice  of  0  •  From 
the  various  functions  examined,  a  hyperbolic  form  was 
selected: 


to’ 


f,  i  -  (ujl) _ y 

;  i+  (tuV/s  / 


(I.2.b*5) 


in  which  the  parameter  A  allows  the  choice  of  the  total 
twist  angle* 

It  should  be  noted  that  the  numerical  work  for  the 
infinitesimal  blade  spacing  is  exploratory  at  best  and  the 
final  design  of  axial  flow  turbomachines  will  be  the  result 
of  the  analysis  involving  a  finite  number  of  blades* 
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One  more  remark  ie  necessary.  The  two-dimensional  vor- 

I 

tex  flow  may  be  obtained  as  a  special  case  of  the  presented 
three-dimensional  analysis.  Then,  v.  and  Z(  are  zero.  It 
was  not  thought  necessary  to  include  numerical  work  for  the 
two-dimensional  analysis,  since  the  choice  of  ur  as  a  con¬ 
stant  would  make  the  flow  a  rather  restricted  one.  However, 
the  solution  of  the  vortex  flow  problem  with  -or  =  -or  (r,x^ 
may  be  obtained  in  a  similar  manner. 

From  the  condition  on  the  selected  shroud  shape,  i.e., 

bxs  =  const.  =  Ms  (I.2.b.6) 

the  shroud  radius  may  be  obtained  as 

rs  -  y\£,  Msz  (I.2.b.7) 

Hence,  £  ,  (x'i  takes  the  form 

Z,  =  'UTo  Ms  (rS‘  *•  Mzz)  (I.2.b.H) 

From  eq.  (I.2.a.l3) 

XL  =  -ur0  Ms  (rs^  i-Ms7.)/r  (1.2. b.9) 

Eq.  (I.2.a.26)  yields  for  the  r  -  Z  variation  of  the 
streamlines 

ra  =  +>  4  Ms  +  M*  z  *■  (I.2.b.l0) 
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? 


From  eq.  (I*2*b*5),  th«  tangential  component  of  the 
▼eloeity  is 

lT=  f/-  ('■£■)}  /{t  +  (I.2.b.ll) 

and  eq*  (I*2*a*26)  yields  for  the  9  -  z.  variation  of  the 
streamlines 

9*  (  wf)  +0  -z}  (I*2.b*12) 

At  the  exit 

Wi)»  (*$r){(e+\)j^[(-k)+\]  -1}  (i.2.b.i3) 

and  this  represents  the  total  twist  of  the  "blade"  (Stream¬ 
line  surface)*  The  arbitrary  constant  fi  allows  the  selec¬ 
tion  of  the  total  twist* 

The  expression  for  the  pressure  distribution  is  obtained 
by  substituting  (I.2*b*l)  into  (I*2.a»24)»  This  gives 

P*  Pi  Wa>4ra-  /  (£0  dz  (I.2.b*14) 

*0 

For  the  case  of  no  pre-rotation  of  the  flow  entering 
the  rotating  stage,  vt<  *  or;  •  IT  the  total  energy 

is  taken  as  constant  for  all  entering  streamlines,  then 

P-  •=  -  (- ur0x+  (I*2,b*15) 

where  denotes  the  ratio  of  the  total  pressure  at 

entrance  to  the  density*  Use  of  (I.2*b*15)  in  (I*2*b*14) 
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yields 


P  *  f*L+  (oi^-v3)/^  (r -*)dz  (I.2.b.l6) 

The  pressure-density  ratio  distribution  at  the  entrance 
is 

Pi  *  ^  -  ¥*"  (I.2.b.l7) 

The  power  imput,  torque,  mass  rate  of  flow,  useful  power 
output  and  the  average  pressure  rise  through  the  stage  may 
be  computed  easily.  However,  such  calculations  are  not  in¬ 
cluded  in  this  example,  since  they  are,  at  best,  approxima¬ 
tions.  They  are  given  in  Ref.  26. 

Table  I  contains  the  numerical  values  used  in  the 
illustrative  example.  Tables  II,  III,  IV,  and  V  give  the 
numerical  results  of  the  flow  parameters  for  the  numerical 
parameters  selected  in  Table  I.  Table  II  gives  the  twist 
<p  computed  for  various  values  of  z  (axial  dimension)  from 
^  ■  0  to  This  table  is  good  for  all  streamlines. 

That  is  to  say,  the  p  -  z  variation  is  independent  of  r 
due  to  the  selection  of  v/r  as  a  function  of  z.  alone. 
Tables  Ilia,  Illb,  and  IIIc  give  the  variation  of  r  ,  xl  , 

V' t  V,  and  P  with  z  along  the  hub,  a  medium,  and 
shroud  streamlines  respectively.  Table  IV  lists  the  varia¬ 
tion  of  static  and  total  pressures  at  exit  for  various 
values  of  •  Table  V  lists  the  variation  of  static 
pressure  at  inlet  for  various  values  of  r;  • 
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The  r  -  2.  variation  of  streamlines  are  shown  in 
Fig.  3»  and  their  -  z  variation  in  Fig.  4*  The  radius 
variation  with  Z  is  presented  in  Fig.  5  in  non-dimensional 
form.  Figs.  6  and  7  show  the  variation  along  z  of  the 
velocity  components  and  of  the  resultant  velocity.  The 
pressure  ratio  distribution  is  presented  in  Fig.  8.  Fig.  9 
is  a  schematic  drawing  of  the  resulting  streamline  surface 
("blade"). 

1.2. c.  ANALYSIS  OF  LAMINAR  OR  TURBULENT  INCOMPRESSIBLE  OR 
COMPRESSIBLE  FLOWS'  WITH  VARIABLE  AXIAL  VELOCITY 
COMPONENT. 

In  the  analysis  that  follows,  the  functions  (1.1.16) 
were  chosen  as 

RsR(r)  and  Z  =  Z  (l)  (1.2.  c.l) 

Eq.  (1.1.12)  multiplied  by  'A  V  yields  as  the  equation  of 
the  r  - z  streamline  trace  (1.1.13) 

R-’Z,  4-  cr  (  <T  »  Ci  »  a  constant)  (I.2.C.2) 

From  (1.1.16)  and  (I.2.C.1),  the  following  relationship  is 
obtained 

2L.  -  7/ 

%r  (I.2.C.3) 

where  prime  denotes  differentiation  (with  respect  to  r  and 
Z  ,  respectively). 
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With  (I.2.C.3),  eq.  (1.1.23)  becomes 

Ztf  fp  rvO  -  (eciO  (I.2.C.4) 

which  again  is  a  Lagrange  differential  equation  with  the 
subsidiary  equations 
sLS,  -  2^3  -  d.(Prv l) 

x7  ^  (I.2.C.5) 

A  general  solution  of  (I.2.C.4)  Bay  be  written  as 

X  (*#-  >  I-Z)  =  O  (I.2.C.6) 

where  X  is  an  arbitrary  function. 

From  (I.2.C.6),  ul  may  be  expressed  in  the  form 

2/ 

u  =  Tr7  X  (R'Z)  (I.2.C.7) 

where  X  is  the  inverse  function  of  X  and  is  determined 
from  the  initial  conditions. 

Eq.  (I.2.C.3)  yields 

T*r  =  $K  XU-Z)  (I.Z.c.4) 

and  'K  V  may  be  found  as 

■fW  =  -TT-X  (R-Z1)  (I.2.C.9) 

In  laminar  and  turbulent  flows,  the  relationship 
between  and  V  is  known 
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H-T  F(v) 


(I.2.C.10) 


Then,  if 

fW  F(v) 

the  expression  for  the  resultant  velocity  will  be 


I.2.C.11) 


X(k-z) 


^  • 


(I.2.C.12) 


where  &  is  a  function  defined  by  (I.2.C.11). 

For  the  special  case  of  incompressible  flow,  the  density 
is  constant,  (°  ■  ■  const.,  and  the  pressure  is  found 

from  eqs.  (1.6),  (1.21)  and  (I.2.C.12)  as 


TP 


=  X(ff-z))  +  iioV>Zj  +0a.((To)  (I.2.C.13) 


along  the  6  -<y-L  streamline  trace. 

For  compressible  flows,  if  adiabatic  change  of  state  is 


assumed  through  the  stage,  or 

A»(*)r 


(I.2.C.14) 


and  the  density  variation  along  any  streamline  trace  may  be 
found  as 


(I.2.C.15) 
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The  pressure,  resultant  velocity  and  Its  components  may 
then  be  found  from  the  given  relationships. 

I.2.d.  NUMERICAL  EXAMPLE  FOR  LAMINAR  COMPRESSIBLE  ADIABATIC 
FLOWS  WITH  VARIABLE  AXIAL  VELOCITY  COMPONENT. 

In  this  example,  the  following  numerical  values  were 

used 

Length  of  stage  ^  ~  0.25  ft 

Constant  angular  velocity  cj  *  1570.8  rad/sec 

Velocity  components  along  leading  edge  of  blade 
Radial  velocity  *Uo»  10  ft/sec 

Axial  velocity  it)©  ■  450  ft/sec 

Tangential  velocity  -  1250  >Tr  ft/sec 

Density  0.002378  slugs/ft^ 

Pressure  2116.8  lb/ft2 

Kinematic  viscosity  *V©»  yr-  *  1.6  x  lO"^  ft2/sec 

Number  of  blades  48 

The  equation  of  the  leading  edge  was  chosen  as 

r0  =  Z  +-  CL  and  r  £  b  (1.2.d.l) 

with  a.  ■  0.25  ft  and  b  *  0.33  ft. 

Fig.  10  shows  the  r  -  z  projection  of  typical  stream¬ 
lines.  Again,  as  in  the  example  of  Chapter  2.b,  the  para¬ 
bolic  distribution  of  the  resultant  velocity  is  used 
(Ref.  25 ).,  see  Fig.  2.  Then,  F(v)  in  e<l*  (1*2. c. 10) 
becomes 
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F(V)*  Cj-Ta. 

with 

C3=  -/2^0 

The  functions  in  (I.2.C.1)  will  be 


(I.2.d.2) 


(I.2.d.3) 


(I.2.d.4) 


(1.2. d.5) 


Fig.  11  shows  the  streamlines  passing  through  r«a_. 
(hub),  v  -  ^3^  (medium)  and  r  *  b  (shroud). 

Fig.  12  presents  the  resultant  velocity,  V  »  variation 
through  the  stage  and  Fig.  13  the  velocity  components.  The 
static  pressure  variation  is  given  in  Fig.  14  and  the  densi¬ 
ty  variation  in  Fig.  15*  From  Fig.  14  it  may  be  seen  that 
the  pressure  rise  across  the  stage  is  approximately  11%. 
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1.2*®.  NUMERICAL  EXAMPLE  FOR  TURBULENT  COMPRESSIBLE  FLOW 
WITH  VARIABLE  AXIAL  VELOCITY  COMPONENT. 

The  initial  values  for  the  turbulent  case  are  the  same 

as  those  for  the  laminar  one  (Chapter  2.d),  except  that 


ir0  ■  674  ft/sec  (I.2.e.l) 


Von  Karmanfs  one-seventh  power  law^  is  used  for  the 
profile  of  the  resultant  velocity  between  blades ,  and 
Schlichting’s  expression  for  the  average  shear  stress-'  • 
These  give 


F  (V)  = 


>%r 


(I.2.e.2) 


with 

C.u-  O.07Z.  (£)^~  (I.2.e.3) 

The  average  value  of  the  resultant  velocity  is 

(X.2.e.4) 

The  results  for  this  example  are  presented  in  Figs.  16- 
20,  showing  the  blade  shape,  resultant  velocity  variation, 
Telocity  components,  static  pressure  and  density  through  the 
stage,  respectively. 
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II.  INCOMPRESSIBLE  VISCOUS  FLOWS  THROUGH  TURBOMACHINE 
STAGES  WITH  LARGE  NUMBER  OF  BLADES.  FINITE  BLADE  SPACING 


II. '1.  METHOD  OF  SOLUTION.  INCOMPRESSIBLE  FLOWS. 


In  this  analysis  the  Lorenz  force-field  is  eliminated 
and  the  three-dimensional,  incompressible  continuity  and 
Navier-Stokes  .equations  are  used.  The  equations  are  used 
in  their  general  form  taking  the  rotation  of  the  stage  into 
account.  However,  the  same  equations  and  analysis  are  valid 
for  stationary  stages  except  that  the  angular  velocity  is 
zero  for  non-rotating  stages. 

The  following  basic  equations  are  used  in  cylindrical 
coordinates: 


The  continuity  equation 

lr)7r  Ir/if  0 

The  Navier-Stokes  equations 

D?  +  +  *r  ^  + 

s.  du.  Cv  -  cor) 

+*  rx  T?"  ”  ~  7* 


(II. 1.1) 


(II.1.2) 
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where 


T  ,  9 
U,  V 

P  «“£■ 

V 

IL  ~  A 

Dt  “  UF 


,  z  are  the  cylindrical  coordinates 

,  ur  are  the  velocity  components  in  the 

y '  ,  <p  *  x  directions,  respec¬ 
tively 

is  the  static  pressure  over  the 
density 

is  the  kinematic  viscosity 
is  the  viscosity 

T  (or+tDr)  +  ixT  &x. 


+  T^  bV* 


The  streamline  equations  for  this  case  are 
_  rd.9  _  d~z- 

UL  “  iT  zxT 


(II. 1.3) 


Infinite  series  in  powers  of  the  tangential  coordinate, 
‘P  ,  are  introduced  for  the  three  velocity  components,  ul  , 
ir  and  nr  ,  and  for  the  pressure  variable,  P  •  The 
coefficients  of  these  series  are  functions  of  the  radial  and 
axial  coordinates,  r  and  7L  •  These  series  are  substituted 
into  eqs<  (II. 1.1)  and  (II.1.2)  and  recurrence  formulae  are 
determined  for  the  coefficients  of  the  series.  Conditions 
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on  the  flow  variable  distribution  along  a  surface  between 
adjacent  blades  as  well  as  boundary  conditions  are  used  for 
the  evaluation  of  the  unknown  coefficients  in  terms  of  which 
the  remainder  of  the  coefficients  of  the  infinite  series  are 
expressed. 

The  number  of  blades  is  chosen  and  a  corresponding  num¬ 
ber  of  streamline  surfaces  are  "frozen".  These  frozen 
streamline  surfaces  are  those  which  have  been  determined  in 
the  analysis  of  the 'turbomachine  stage  with  infinitesimally 
spaced  blades  (Chapter  I.2.b).  The  distribution  of  the 
flow  variables  along  these  "frozen"  streamline  surfaces  are 
assumed  to  be  the  same  as  those  calculated  in  the  "infinites¬ 
imal"  analysis  and  are  denoted  by  subscript  zero.  If  the 
first  coefficients  in  the  four  series  are  taken  to  be  equal 
to  the  corresponding  values  of  the  flow  variables  along  the 
frozen  streamline  surface,  the  flow  variables  will  vary  from 
these  values  toward  the  actual  blade  surfaces  of  the  "finite" 
case.  Therefore,  the  tangential  variable  is  measured  from 
the  "frozen"  streamline  surface.  The  introduction  of  these 
"frozen"  streamline  surfaces  and  the  prescribed  distribution 
of  the  flow  variables  along  them  will  serve  to  ensure  a  more 
rapid  convergence  of  the  series  solution. 

The  boundary  conditions  used  for  the  evaluation  of  the 
remaining  unknown  coefficients  in  the  series  are  the  non¬ 
slip  conditions  along  the  surfaces  of  the  blades  and  the 
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closure  of  the  streamlines,  which  form  the  blades,  at  the 
leading  and  trailing  edges* 

The  blades  of  modern  turbomachines  are  thin  and,  there¬ 
fore  ,  it  is  thought  to  be  suf f icient  to  enforce  the  non-slip 
conditions  along  the  chord  surfaces  of  the  blades  instead  of 
at  the  surfaces  themselves.  This  relaxation  of  the  non-slip 
conditions  greatly  simplifies  the  analysis  and  the  numerical 
computations,  since  the  enforcement  of  the  non-slip  condi¬ 
tions  along  the  as  yet  to  be  determined  blade  surfaces  would 
make  the  boundary  conditions  highly  non-linear.  The  error 
introduced  by  this  relaxation  is  thought  to  be  negligible. 

The  closure  of  the  blades  at  the  leading  edge  is  ensured 
by  chosing  the  streamline  surfaces  forming  both  the  Upper  and 
the  lower  sides  of  the  blades  so  as  to  pass  through  the  lead¬ 
ing  edge.  This  is  possible  by  the  choice  of  two  arbitrary- 
constants  which  arise  from  integration  of  the  streamline 
equations  (II.1.3).  The  closure  of  the  blade  surfaces  at  the 
trailing  edge,  however,  cannot  be  achieved  in  general  form, 
but  has  to  be  done  for  each  design  case  by  judicious  choice 
of  two  arbitrary  functions  representing  two  of  the  three 
coefficients  of  the  first  power  terms  of  the  tangential 
variable  in  the  series  for  lx  ,  v*  and  ~ux  . 

The  series  for  the  velocity  components  are  taken  in  the 

form 
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(II. 1.4) 


a6=  0-^.S au5n] 

■n.*/ 

where 


t)/ »  <p  -  <pfl  (-r,x)  is  the  tangential  variable 

measured  from  the  "frozen" 
streamline  surface 


is  the  tangential  coordinate 
of  the  "frozen"  streamline 
surface 


is  the  tangential  coordinate 
measured  from  <pa  to  the  blade 
chord  above  <p0 


is  the  tangential  coordinate 
measured  from  %  to  the  blade 
chord  below  <P0 


denotes  any  of  the  velocity 
components  xl  ,  xr  ,  or  ixr 


denotes  the  corresponding 
"infinitesimal"  velocity  com¬ 
ponents  xl  $  t  TTo  ,  or  -ur0 


?o  (r,z) 

-y/u 

OL 

a-o 
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? 


Qu-n  denotes  the  coefficients 

in  the  series  for  the  velocity 
components,  ir^  ,  or  -ur^_ 


The  series  for  the  pressure  variable  is  chosen  as 
Te  =  [j  ,  P-nS71]  (II. 1.5) 

where 


P0  is  the  pressure  variable  along  the  "frozen" 
streamline  surface 


is  the  coefficient  in  the  series  for  the 
pressure  variable 

V 

For  Y*  *  0,  5*0  and  the  expressions  for  the 

velocity  components  and  the  pressure  variable,  (II. 1.4)  and 
(II. 1.5),  satisfy  the  conditions  that  along  the  "frozen" 

-4  streamline  surface  CL*  cl0  and  P-  P0  .  Since  V  ■ 

and  V/*VW’  denote  the  two  chord  surfaces  of  the  blades 
adjacent  to  the  "frozen"  streamline  surface,  it  can  be  seen 
from  the  general  expression  of  the  velocity  components 
(II. 1.4)  that  for  these  values  of  Y'  bhe  velocity  components 
,  will  be  zero.  Therefore,  the  non-slip  conditions  along  the 

i  (  chord  surfaces  will  be  satisfied  by  the  expressions  for  the 

velocity  components. 

Fig.  21  shows  the  notation  used  in  the  9  -  Z.  plane. 

i 

From  this  figure  it  can  be  seen  that 
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(II. 1.6) 


where 

is  the  angular  distance  between 
the  mean  .line  and  the  "frozen" 
streamline  surface;  the  mean 
line  lies  half  way  between  two 
adjacent  chord  surfaces 

Since  the  blade  spacing 

<| +  (II. 1.7) 

the  choice  of  ^4  and  4  determines  *Y^  ,  and  £  • 

Eq.  (II. 1.4)  written  for  the  three  velocity  components  be¬ 
comes 

Ue"*  [/-SfcS-S'tJC  I  +*  f  lO71] 

"ns/ 

-  [/ -sa]L  i  +  1  v^sn] 

•n.«l 

5  -St)C  I.  +  -bf-n.  S U] 

By  substituting  eqs.  (II. 1.5)  and  (II. 1.8)  into 
eqs.  (II. 1.1)  and  (II. 1.2)  collecting  terms  of  equal  powers 
of  S  and  equating  their  coefficients  to  zero  one  obtains 
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four  simultaneous  algebraic  equations  for  each  set  of  four 
coefficients.  The  solution  of  these  equations  yields  the 
recurrence  equations  for  the  coefficients  of  the  series  in 
terms  of  U0  ,  V'o  ,  \Uo  ,  P0  ,  \l,  and  either  ift  or 
"UT)  •  The  only  two  arbitrary  functions  of  r  and  z  amongst 
these  are  ix,  and  either  y,  or  iy,  •  These  two  arbitrary 
functions  may  be  used  to  ensure  closure  of  the  blade  surfaces 
at  the  trailing  edge. 

Before  presenting  the  recurrence  equations,  a  brief 
discussion  of  Y  and  its  partial  derivatives  should  be 
given.  After  eq.  (II. 1.4),  Y  is  defined  as 

V*  9  -V.*)  {ii. 1.9) 

It  can  be  seen  from  eq.  (II. 1.9)  that  &V/69  equals 
unity,  but  both  r  and  are  functions  of 

T  and  Z  •  If,  however,  the  "frozen"  streamline  surface 
is  made  up  of  radial  lines,  i.e.,  if  90  ~  only, 

then  is  unity,  is  zero  and  the 

is  ~V'0/tvJ'o  »  as  can  be  calculated  from  the  equations  of 
the  streamlines,  eq.  (II.1.3),  for  the  "infinitesimal"  case. 
In  this  somewhat  restricted  case,  the  derivatives  of  the 
tangential  variable  5  simplify  and  the  recurrence  equa¬ 
tions  for  the  coefficients  of  the  series  become  much  more 
manageable.  After  eq.  (II. 1.4),  $  is  defined  as 
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(II.1.10) 


(II.l.ll) 


It  should  be  noted,  however,  that  the  restriction  on  the 
"frozen"  streamline  surface  that  it  be  made  up  of  radial 
lines  does  not  necessarily  restrict  the  blade  surfaces  of  the 
"finite"  case  to  radial  ones* 

In  the  following,  the  recurrence  equations  are  presented 
for  the  general  case,  i.e.  for  the  case  where  9  =  9 ) 

The  recurrence  equations  give  general  expressions  for  the 
coefficients  from  ul8  ,  Vg  ,  ~UTt  and  P7  on  (including 
these  coefficients). 

The  recurrence  equations  are 

=  (II. 1.12) 

Vun ♦  \  R*-,  =  -  *£„  dl.1.13) 


3  =  yj/fpryK. 

and  its  partial  derivatives  are 
4r  =  0 

is  _  1 

69  ~ 


b_S  _ 
b  x  ~ 


~VL 
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SCnwn  %  Pn-i  =  "  t-n  (II. 1.15) 

/ 

In  the  above  four  equations,  the  expressions  for  u.-*.  , 
"W  ,  and  P-n  are  determined  from  the  following 

determinants.  The  fundamental  determinant  is  denoted  by 
(F.D.) 

“A-n  X  O 

0  0  *0* 

0  r6-n  0  X 

0  0  £0„ 
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(F.  D.)irn  ~ 


(FO.^?iru  = 


^o.)a,.r- 


% 

o 

V 

*E 

0 

*Dn 

0 

% 

0 

*>n 

0 

%* 

*-77 

% 

"On 

*A 

•^-n 

O 

A-rv 

0. 

c'n 

fa* 

0 

K 

*E 

fa* 

0 

0 

rE 

c>v. 

fan 

<*c 

v-n 

Ot- 

c_n. 

V 

0 

0 

'*£n 

0 

yB. 

0 

0 

0 
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The  series  for  "lx  ,  ir  and  ur  allow  the  numerical 


integration  of  the  streamline  equations 

cLt  _  rot  9  _  cLz 
il  “  ir  ”  zlT 


(II. 1.21) 


For  the  determination  of  the  r  and  projections  as 
functions  of  i  ,  eq.  (II. 1.21)  will  be  integrated  by  Taylor 
series  expansions  about  consecutive  points  of  z.  along  the 
streamlines. 

The  expressions  used  in  eqs.  (II. 1.12)  to  (II. 1.15)  are 
quite  lengthy  and  are  not  included  here.  They  are  presented 
in  Ref.  28, 


11. 2.  NUMERICAL  WORK 

11. 2.  a.  PRELIMINARY  REMARKS 

In  the  general  analysis  of  laminar,  incompressible  flows 
through  turbomachines  with  finite  blade  spacing  and  finite 
blade  thickness  (Section  II. 1),  the  velocities  and  the  pressure 
density  ratio  are  expanded  in  infinite  series  in  the  tangential 
variable  9  about  points  9=  9©  ,  where  {90('r,x')l  is  the 

set  of  points  on  the  "frozen"  streamline  surface  (see  Fig.  21). 
The  recurrence  equations  determine  the  Coefficients  of  the 
series  in  terms  of  ix.e('r)z),  ■v/0(r,z)  ,  -ur0('c>z),  ix,(r,z\ 

and  (or  itr,(r,i)  ,  since  iry-yrl  by  the  first  recur¬ 

rence  equation  which  is  derived  from  the  continuity  equation). 
The  functions  il0  ,  Ve  »  "Wo  and  Pe  are  known  from  the 
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analysis  of  incompressible  flows  through  turbomachines  with 
infinitesimally  spaced  blades  (Section  I.1J.  The  "frozen 
streamline"  surface,  <p  =  <p0  (  r,0  ,  and  the  approximate 

blade  shape  are  also  known. 

There  remain  to  be  chosen  only  the  functions  u.,  (v.-zS) 
and  ir;  (r.-O  .  These  functions  are  thought  to  be  deter¬ 

mined  by  the  conditions  of  closure  of  the  blade  surface  and 
of  convergence  of  the  infinite  series. 

II. 2. b.  THE  BLADE  SHAPE  -  CLOSURE. 

The  numerical  example  of  Section  I. 2. a  with  16  blades 
has  as  solution 

“U0  =  ~ur0  Ms  (rs.  -H  Msz)t 

Vo  -  T  (O  \  I  +  jti  ) 


WQ  —  constant 


+  J" 


4$ 
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-a  m[-  <*> 


(II.2.b.l) 


_  rf  -  Maz(r3i  tMsz)  -  MszUrsi+M^-r*]* 

C  -  r4  -  MS7.  (ru+  Msx) ■>■  Msz.[(riC-h Ms7.')-rJJz 

For  this  solution,  the  velocity  profile  was  chosen  as  a 
parabola  with  zero  eccentricity.  Having  chosen  this  symmet¬ 
ric  profile  means  that  the  "frozen  streamline",  f  -  <P0(r,-x')  , 
may  be  fixed  in  the  middle  of  the  blade  channel,  i.e., 

'vVu='ifJ-<-  ,  so  that  HH-O  in  eqs.  (11.1*4) • 

The  problem  may  be  stated  as  the  determination  of  the 
as  yet  unknown  functions  n/r.x)  and  zr,  (r,z)  such  that 
the  infinite  series  are  convergent  in  the  region  between  the 
blades  and  that  the  blade  of  finite  thickness,  determined  by 
the  velocity  series,  closes  at  the  leading  and  trailing  edges. 
The  form  of  the  series  assures  that  the  components  UL  ,  IT  , 
IlT  and  P  will  reduce  to  ilq  ,  ir0  ,  \tr0  and  P©  on  the 
"frozen  streamline"  surface,  •  The  requlre- 
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ment  that  the  velocity  is  a  maximum  on  this  surface  may  be 
met  by  enforcing 


(II.2.b.2) 

and 

wL=o 

Which,  upon  substituting  the  velocity  series,  eq.  (II. 1.4), 
clearly  implies  that  the  coefficients  =  xt;  -  O 

and  also  A.  *  0  for  all  values  of  -r  and  z.  .  It  is  not 
immediately  clear,  however,  that  this  choice  of  the  arbitrary 
functions  Tjl,  and  xr,  will  lead  to  closure  of  the  blade 
surface  and  convergence  of  the  series. 

There  are  several  ways  to  treat  the  condition  of  closure. 
One  of  the  most  natural  of  these  is  to  specify  the  leading 
edge  of  the  blade,  determine  the  functions  xc,  and  xr, 
from  other  considerations,  as  for  example,  the  maximal  con¬ 
dition  of  the  preceding  paragraph,  and  allow  the  resulting 
blade  to  close  where  it  will,  if  at  all.  If  the  resulting 
blade  does  not  close,  it  is  necessary  to  reconsider  the 
functions  xc,  and  v;  ,  alter  them  in  some  way,  rework  the 

numerical  example  and  examine  the  resulting  blade  for  closure. 
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This  procedure  must  be  repeated  until  closure  is  attained* 

I 

However,  to  simplify  the  design  of  the  stator  stage,  it  is 
desirable  to  formulate  the  closure  condition  so  that  both 
the  leading  and  trailing  edges  of  the  blade  of  the  rotor 
are  specified*  This  may  be  done  as  follows: 

The  upper  blade  surface,  9-  9^ (v/i)  ,  is  generated 

as  the  solution  of  the  equation 


ix  ^5*  +  -  T"  (11*2. b*3) 

such  that  meets  the  conditions  that  where  X  *  0, 

^u.  (t,  o)  =  %  ( r\o)  -t- "9^  and  where  , 

^-u.  «,  Similarly,  the  lower  blade 

surface,  <p  -  q?  (<^7.) ,  is  generated  as  the  solution  of 


(II.2.b*4) 


such  that  where  1=0,  %(r,o)~  90(r> o')  -  and  where 

-  9,  (r,^  "  ‘9i  It  should 

be  noted  that  the  subsidiary  equations  of  (II.2*b*3)  and 
(II*2*b*4)  are  the  streamline  equations,  i.e. 


Sir  -  -  rT  dL9 

■u.  ~  xo=  ir 


(II.2.b*5) 


Here,  the  leading  and  trailing  edges  are  specified  to  be 
radial  lines  with  the  same  "twist"  as  the  blade  of  the  in¬ 
finitesimal  case*  This  statement  of  the  problem  may  be 
written  in  terms  of  line  integrals  along  the  characteristics 
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(streamlines)  of  eqs.  (II.2.b.3)  and  (II.2.b.4).  The  appro¬ 
priate  integral  representation  is 


(II.2.b.6) 


where  V,  Is  a  streamline  passing  through  any  initial  point 
1  =  0  »  T-v,i  ,  on  the  initial  line  9uJ<\o'>=9e(r>o')+'H)U' 

(or  =  90  (r.o')  -  ^  ),  terminating  in  the  plane 

3L  -  JL  •  /^o  is  the  total  angle  of  twist  of  the  "infini¬ 


tesimal”  blade. 

From  1. 2. a  and  I.2.b, 


A- 


-jrn 

rur0  cLz. 


(II.2.b.7) 


where  is  a  streamline  passing  through  any  initial  point 

4  =  0  ,  r  =  r '  and  terminating  in  the  plane  X-  =  . 

'Uo 

Since,  from  eqs.  (II.2.b.l),  the  ratio  r^r0  is  independent 
of  r  and  <P  ,  the  "infinitesimal"  blade  is  a  surface  gener¬ 
ated  by  radial  lines  and  /3e  is  a  constant  (independent  of 

irQ 

r  ).  It  is  also  to  be  noticed  that,  since  the  ratio  -rure 
is  independent  of  r  and  9  ,  the  value  A>  is  independent 
of  the  path  of  integration  between  the  planes  X-O  and 

7--.Q  .  So  the  streamline,  ,  may  be  replaced  by  any 

path  extending  from  the  entrance  plane,  7.-0  ,  to  the 

exit  plane  Z-  JL  •  In  particular,  it  may  be  replaced  by 
the  curve  V7,  • 
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Combining  (II.2.b.6)  and  (II.2.b.7),  one  sees  that 


TUT  cLt.  =A 


-  fas. 

“  I  ThTo  dLZ 


f  --o 

Jri 


(II.2.b.S) 


is  a  statement  of  the  condition  of  closure. 


Substitution  from  (II. 1.4)  into  (II.2.b.7)  yields 


/  +2  uViS*  1  5  wi  r-  ~  O 


(II.2.b.9) 


This  condition  is  difficult  to  manipulate  in  its  present 
form.  If,  however,  it  is  desired  to  enforce  the  closure 
condition  to  only  the  second  order  approximation  of  the 
velocity  series,  so  that  (II.2.b.9)  may  be  written 


f  { l+ifiS+T&S.*  .  .'I 
/  ( I  s  I  S 


arc,  dx  _ 
urc  -r  ~  O 


(II. 2. b. 10) 


it  becomes  possible  to  find  elementary  conditions  on  tx, 
and  ir,  (or  ur,  )  such  that  closure  occurs  to  this  approxi¬ 
mation. 

Clearing  fractions  in  the  integrand  yields 
I  i  *  ur,  5  +  5^  r  ~0  (II, 2. b. 11) 


From  the  first  continuity  recurrence  equation  and  the 
functions  in  (II.2.b.l)  it  is  seen  that 
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\K  (rX)  =  "ur,  (• r,z.) 


(II. 2. b. 12) 


so  that 


(II. 2. b. 13) 


/jfa-frtOs*-  V5 

7  /  +  xir,5  We  r  u 

From  the  recurrence  equations  of  Section  II. 1, 

”‘84(’Si'Di-7DJ'0  (II. 2. b. 14) 


and 


...  ,  *cJbA  ♦"Bifefo-EjpQ 


(II. 2. b. 15 ) 


if  °%a-o 


so  that  'u^=iir4< 

Using  the  special  forms  (II.2.b.l),  °tEi-0  reduces 
to 


ii-o  *  xw-o  ^  =  o 


(II.2.b.l6j 


The  integral  in  (II. 2. b. 13)  is  zero  if  'V'^-'ur^  ,  and  the 

closure  condition  may  be  met  by  choosing  and  Vi  (-  -urt) 

such  that  (II.2*b.l6)  is  satisfied. 

As  was  pointed  out  in  the  discussion  of  eq.  (II.2.b.2), 
the  condition  that  a  maximum  of  the  velocity  occurs  on  the 
surface  <P=  %  (V,x)  ,  may  be  met  by  choosing  the  functions 

i  and  v*,  zero.  Such  functions  trivially  satisfy 
eq.  (II.2.b.l6),  (since  ir,  =  ur(  )  and,  therefore,  were 
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tried  for  the  numerical  example.  Unfortunately,  the  result¬ 
ing  velocity  series  appear  to  be  divergent.  Several  possible 
reasons  for  this  will  be  discussed  in  Section  II.2.C.  But 
now,  it  appears  proper  to  discuss  eqs.  (II. 2. b. 13)  and 
(II.2.b.l6)  in  slightly  greater  detail,  and  show  that  eq. 

(II. 2. b. 16)  leads  to  a  thin  airfoil  shape. 

If  the  integrand  of  eq.  (II. 2. b. 13)  is  identically  zero 
(  vz  ■=  zlT x.  )#  then,  to  the  second  order  approximation,  the 
"blade"  must  close  everywhere,  i.e.,  the  integral  (II. 2. b. 13) 
becomes  a  function  identically  zero  in  f  ,  z  and  S  .  This 
means  that  the  functions  resulting  from  the  integration  of 
(II.2.b.3)  and  (II.2.b.4),  <PU  and  9*  (where  the  series 
coefficients  in  these  equations  have  been  terminated  at  the 
term  in  5A  )  must  differ  only  by  the  angular  distance  be¬ 
tween  successive  "frozen"  stream  surfaces,  i.e. 

=  Y*-yA  (n.2.b.i7) 

This  condition  corresponds  to  a  blade  of  zero  thickness. 
Further,  the  only  admissible  blade  of  zero  thickness  is  the 
original  "chord  surface",  9  =  90  (V.z')  +‘'f’Xc  or 
9=  9e  ( v , tS)  -  y*  •  For,  suppose  there  were  an  admis¬ 
sible  blade  of  zero  thickness  different  from  this  blade, 
then  such  a  blade  must  have  points  above  or  below  the  orig¬ 
inal  chord  surface  upon  which  the  velocities  are  zero.  And, 
since  the  blade  is  generated  as  a  family  of  streamlines,  the 
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flow  must  somewhere  have  passed  through  the  original  chord 
surface*  But  flow  with  zero  velocity  is  stationary  and  so 
could  not  pass  through  this  surface.  A  contradiction* 

It  is  seen  that  if  the  closure  condition  (II*2*b*13)  is 
enforced  identically  by  making  ^EL^-O  ,  then  any  devia¬ 
tion  of  the  final  blade  shape  from  that  of  the  original  chord 
line  airfoil  must  be  due  to  the  influence  of  terms  in  the 
velocity  series  of  order  5 3  and  higher*  Such  an  airfoil 
must  then  be  a  thin  airfoil* 

II. 2. c*  THE  NUMERICAL  EXAMPLES 

Before  proceeding  to  a  detailed  consideration  of  the 
numerical  examples ,  it  should  be  pointed  out  that  derivatives 
of  preceding  coefficients  up  to  ,u.n«,  „  > 

and  P.„_|  occur  in  the  recurrence  equations  for  'LL-*., 

,  ur-n.  and  Pu  •  The  calculation  of  these  derivatives 
is  a  laborious  process  which  was  carried  out  only  for  the 
specific  examples  considered*  Because  of  the  special  nature 
of  these  rather  lengthy  calculations,  they  are  not  included 
in  this  report,  but  are  placed  on  file  in  the  departmental 
office  in  Urbana  and  are  available  there  to  interested  per¬ 
sons* 

Following  eq.  (II. 2. b. 16)  in  the  discussion  on  closure, 
it  was  mentioned  that  an  example  has  been  worked  for  the 
solution  of  the  case  with  infinitesimal  blade  spacing,  and 
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for  the  choice  of  the  arbitrary  functions,  *u.(  ■=•  v-)  =  ur,  -  Q 

« 

and  also  "4*0.  Also  noted  there  was  the  fact  that  this 
choice  of  the  arbitrary  functions  leads  to  closure  in  the 
second  order  approximation  and  to  a  velocity  profile  between 
the  blades  which  has  a  relative  maximum  (or  a  minimum)  on  the 
frozen  streamline  surface#  To  guarantee  the  relative  maximum, 
the  slope  of  the  tangent  must  be  proceeding  from  positive 
values  to  negative  values  in  the  neighborhood  of  the  points 
on  this  surface#  This  is  the  usual  condition  of  negativeness 
on  the  second  derivative,  which  in  this  case  is  seen  to  be 
satisfied  if 


Q <.  I  (II#2.c#l) 

where  Clx  represents  each  of  tjl^  ,  1 and  'ur^.  » 

More  important  was  the  fact  that  the  first  three  terms 
of  the  velocity  series  were  seen  to  be  rapidly  divergent  in 
the  region  where  5  =■  *  '2S  #  Since  approximation  formulae 

were  used  for  the  term  P0  ,  as  discussed  in  Section  I«2.b, 
involving  the  assumption  that  **  Brt  9  i#e.  the  stagna¬ 
tion  enthalpy  is  constant  throughout  the  stage,  and  since 
the  coefficients  Cl-n.  depend  on  the  derivatives  ft>,r  , 

9C, 2.  >  Pe.rr  »  Pe,r*  ,  Pc.xa.  ,  etc#,  it  is  thought  that 

the  divergence  of  the  series  is  traceable  to  errors  intro¬ 
duced  by  this  approximation#  Another  possibility,  however, 
is  that  choosing  ix, -ir,  -O  might  be  inconsistent  with 
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the  series  expansion  and  the  differential  equations.  To 
investigate  this  possibility,  a  new  function  xt,  was  chosen 
satisfying  two  conditions;  first,  xl,  =■  xl  ((x)  so  that 
eq.  (II.2.b.l6)  is  satisfied,  and  second,  xl,  (T-')  makes 
K-a  Cr.x')  -  O  at  approximately  the  center  of  the  stage 
(  'C  *0.35)*  This  choice  is  accomplished  by  solving  the 
equation  (  r  =0.35,  z.  )  ■  0,  since  it  can  be  seen 

A- 

from  Section  II. 1  that,  if  E a.  »  0,  then  xc^s  o.  A 
solution,  XL,  (z)  ,  for  this  equation  (with  the  approxi¬ 

mations  on  P0)r  retained)  is  displayed  in  Fig.  22. 

With  this  new  choice  of  to.,  ,  a  second  numerical  cal¬ 
culation  was  then  performed.  Here  again  the  first  three 
terms  of  the  velocity  series  were  noticed  to  be  rapidly 
divergent.  Figs.  23-29  illustrate  this  second  case. 

To  get  an  idea  of  how  the  approximation  on  PG  for  the 
second  case  (xl,‘^o')  t  affects  the  values  of  the  deriva¬ 
tives  Po,r  •)  PojJ.  Pe.rr  ^  P0,r7.  etc.,  Po  ,  Po.r  , 

and  P0)X  were  programmed  for  Illiac  computation  utilizing 
the  complete  solution  X.-XCk'.z)  eq.  (II.2.b.l).  A  cursory 
examination  of  the  formula  for  indicates  that  this 

value  will  be  very  much  larger  than  the  approximated  values 
near  the  shroud  and  near  the  exit  of  the  compressor  stage. 

The  exact  calculations  show  that  the  approximated  values  at 
Z.  ■  0.2  of  P0lr  are  off  by  about  20  per  cent,  and  those 
of  Po>z  by  about  50  per  cent.  The  values  of  these  first 
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derivatives  at  smaller  values  of  Z  appear  to  be  less  in 
error.  One  naturally  expects  higher  order  derivatives  to  be 
more  in  error  than  these  first  order  derivatives.  However, 
the  exact  numerical  values  for  the  higher  derivatives,  e.g., 
Po.xz.  »  Po.rr  >  Po.cjr  *  were  not  computed  with  the 
Illiac,  although  the  formulae  for  their  calculation  have  been 
prepared. 

However,  it  is  not  known  whether  convergence  of  the 
series  will  be  uniformly  effected  by  this  refinement  in  cal¬ 
culation.  It  was  thought  that  a  re-examination  of  the  phys¬ 
ical  principles  underlying  the  choice  of  the  arbitrary  func¬ 
tions  U-,  and  ir,  would  be  helpful.  The  gross  results  of 
the  numerical  examples  already  calculated  are  useful  for 
this  re-examination. 

It  is  evident  from  the  example  calculated  for  tr,  »  q 
and  u,5U,(i )  (such  that  (0.35,  *-  )  **  0)  that  where 

is  small,  i.e.,  in  the  neighborhood  of  r  ■  0.35,  the 
other  coefficients  are  small  also  (c.f.  Fig.  2 4).  This  leads 
one  to  expect  that  convergence  of  the  series  may  be  obtained 
by  choosing  -ul(r,z')  and  tt,  0,-0  such  that  ix^-0 
and  such  that  closure  occurs  to  the  second  order.  It  is 
actually  possible  to  choose  ,  =  v,  =  O  at  the  entrance  of  the 
first  stage.  This  means  that  the  axially  symmetric  flow 
leaving  the  entrance  diffuser  is  undistrubed  to  the  first 
order. 
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The  condition  of  closure  is  given  in  (II.2.b.l6)  and  the 

—  A _ 

condition  O  may  be  met  by  taking  E-&-  O  or 

by 


A 


+  C.  or,  +  0  -  O 


(II.2.C.2) 


where  A  ,  B  ,  C.  ,  and  D  are  functions  of  r  and  x 
known  from  the  zero  order  terms.  Dividing  (II.2.C.2)  by  C  , 
differentiating  with  respect  to  x  and  using  (II.2.b.l6),  one 
obtains  a  linear,  parabolic,  second  order  partial  differential 
equation  in  ix., 


(II.2.C.3) 


The  characteristics  of  this  equation  are  lines  parallel  to 
the  X  -axis.  Therefore,  physically,  a  solution  is  expected 
to  exist  reducing  to  arbitrary  values  of  ul,  and  its  normal 
derivative  on  the  initial  value  curve  X  ■  0.  Eq.  (II.2.C.2), 
with  specification  e.g.,  o  on  the  curve  x  -  0, 

leads  to  a  solution  of  (II.2.C.3)  for  'll,  throughout  the  com¬ 
pressor.  Knowledge  of  xl,  and  again  use  of  eq.  (II.2.C.2) 
yields  7 r,  throughout  the  compressor  stage.  Therefore,  it 
is  thought  that  the  above  procedure  in  solving  the  problem 
may  yield  convergent  series  for  the  flow  variables. 
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III.  VISCOUS  COMPRESSIBLE  FLOWS  THROUGH 
TURBOMACHINE  STAGES  WITH  LARGE  NUMBER  OF  BLADES. 

FINITE  BLADE  SPACING. 

III.l.  LAMINAR  FLOWS. 

The  basic  equations  used  are  the  Navier-Stokes  equations 
together  with  the  continuity  equation  and  the  energy  equation 
in  their  most  general  form.  Infinite  series  in  powers  of  the 
tangential  coordinate  (with  the  coefficients  being  functions 
of  the  radial  and  axial  coordinates)  are  introduced  for  the 
dependent  variables.  The  Navier-Stokes  equations  and  the 
energy  equation  are  of  second  order  and,  therefore,  the 
coefficients  of  the  series  are  expressed  in  terms  of  the 
coefficients  of  the  zeroth  and  first  order  terms  of  the 
tangential  coordinate. 

The  analysis  of  axially  symmetric  compressible  flow 
through  a  rotating  turbomachine  stage  with  infinitesimal 
blades.  Section  I,  gives  the  zeroth  order  terms.  The  first 
order  terms  are  specified  from  the  boundary  conditions  which 
are  the  non-slip  condition  on  the  surface  of  the  blade  and 
the  closure  condition. 

The  method  of  solving  the  problem  using  the  analysis 
for  the  infinitesimally  spaced  blades  and  the  boundary 
conditions  is  explained  below. 
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By  using  cylindrical  polar  coordinates,  the  basic 
equation  can  be  written  as? 

The  continuity  equation 

f>Q  +■  D^bt  =  o 


The  Navier-Stokes  equations 
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/>%¥  =  -'H  >vV-  3fc©(l^’) 


(III. 1.4) 


The  energy  equation 


„  DT.  -  -D£  4.  J- 
/^C,  Dt  Dt  f  r 


a[rA(|g). 
d  r 


+  T* 


+z<i&r+(&-z  ♦+$*)*■ 

wh^ms^w] 


The  equation  of  state 

P=R  PT 


where 


are  the  cylindrical 


(III. 1.5) 


(III.1.6) 


coordinates 
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•u  jvsttr 

P 

/° 

T 

r 

X 

R 


are  the  velocity  components  in 
the  r  {  <9  f  z.  directions 
respectively 

is  the  static  pressure 

is  the  density 

is  the  temperature 

is  the  coefficient  of 
viscosity 

is  the  coefficient  of  heat 
conduction 

is  the  gas  constant 


CP 


is  the  specific  heat  at  con¬ 
stant  pressure 


is  the  angular  velocity  of 
the  rotating  stage 


]>t  ~  Tx(chr)  ( ir*  cor)^  +  icr^r 


The  quantities  X  and  X  are  given  by 
V‘  =  C.T  >  X-oiT  (  C  ,  oL  constants)  (III.  1.7) 

and  Cp  is  constant. 


(III. 1.6) 


The  streamline  equations  for  this  case  are 

-Str  -  rd<P  _  d.z. 

ir  ~  -ur 

The  pressure  P  is  eliminated  from  eqs.  (111.1.2  to  III.1.5) 
by  using  eq.  (III. 1.6).  Infinite  series  in  powers  of  the 
tangential  coordinate  (  9  )  are  introduced  for  the  three 
velocity  components  (  'll  #  tr  ,  ~ur  ) ,  for  the  density  J° 
and  for  the  temperature  T  •  The  coefficients  of  these 
series  are  functions  of  the  radial  and  axial  coordinates 
(  T  ,  z  ).  These  series  are  substituted  into  eqs.  (III. 1.1 
to  III. 1.5)  and  recurrence  formulae  are  determined  for  the 
coefficients  of  the  series.  Conditions  on  the  flow  variable 
distributions  along  a  surface  between  adjacent  blades  as  well 
as  the  boundary  conditions  are  used  for  the  evaluation  of  the 
unknown  coefficients  in  terms  of  which  the  remainder  of  the 
coefficients  of  the  infinite  series  are  expressed. 

The  number  of  blades  is  chosen  and  a  corresponding  num¬ 
ber  of  streamline  surfaces  are  "frozen".  These  frozen  stream¬ 
line  surfaces  are  those  which  have  been  determined  in  the 
analysis  of  the  turbomachine  stage  with  infinitesimally 
spaced  blades.  The  distribution  of  the  flow  variables  along 
these  frozen  streamline  surfaces  are  assumed  to  be  the  same 
as  those  calculated  in  the  "infinitesimal"  analysis  and  are 
denoted  by  the  subscript  zero.  If  the  first  coefficients  in 
the  five  series  are  taken  to  be  equal  to  the  corresponding 
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values  of  the  flow  variables  along  the  streamline  surface* 
the  flow  variables  will  vary  from  these  values  toward  the 
actual  blade  surfaces  of  "finite”  spacing*  Therefore,  the 
tangential  variable  is  measured  from  the  frozen  streamline 
surface. 

The  boundary  conditions • used  for  the  evaluation  of  the 
remaining  unknown  coefficients  in  the  series  are  the  non-slip 
conditions  along  the  surfaces  of  the  blades  and  the  closure 
of  the  streamlines  which  form  the  blades  at  the  leading  and 
trailing  edges. 

The  blades  of  modern  turbomachines  are  thin  and,  there¬ 
fore,  it  is  thought  to  be  sufficient  to  enforce  the  non-slip 
conditions  along  the  chord  surfaces  of  the  blades  instead  of 
at  the  surfaces  themselves.  This  relaxation  of  the  non-slip 
conditions  greatly  simplifies  the  analysis  and  the  numerical 
computation  since  the  enforcement  of  the  non-slip  conditions 
along  the  as  yet  to  be  determined  blade  surfaces  would  make 
the  boundary  conditions  highly  non-linear.  This  same 
“relaxing"  of  the  boundary  conditions  was  used  in  previous 
Sections. 

The  series  for  the  velocity  components  are  taken  in 
the  form 

=  1,  ^S'] 
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(III. 1.9) 


v%  =  [l  '  KS  -SA][  / 1  £  w  s’1] 
Idj  ~[|-KS-s][l 


and,  the  series  for  the  temperature  and  the  density  are 
chosen  as 


"  11 

*  'h 

C A  P 

*  (Jt 

1 

(III.1.10) 

where 

<V-rul 

r\  -  y.u. 

5  ” 

is  the  tangential  variable 
measured  from  the  "frozen" 
streamline  surface 

90(r,z) 

is  the  tangential  coordinate 
of  the  "frozen"  streamline 
surface 
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y  is  the  tangential  coordinate 

measured  from  9s  to  the 
blade  chord  above  ?0 

is  the  tangential  coordinate 
measured  from  *P0  to  the 
blade  chord  below  To 

For  'f'  *0,  5*0  and  the  expressions  for  the  veloci¬ 
ty  components,  temperature  and  density  satisfy  the  conditions 
that  along  the  frozen  streamline  surface  "U.  =  u.  0  , 

ir  -  iTo  ,  etc.  Since  9  =  “  and  9  -  T  ^  denote  the 
two  chord  surfaces  of  the  blades  adjacent  to  the  frozen 
streamline  surface,  frofa  eq.  (III. 1.9),  it  can  be  seen  that 
for  these  values  of  T  the  velocity  components  will  be  zero. 
Therefore,  the  non-slip  conditions  along  the  chord  surfaces 
will  be  satisfied  by  the  expressions  for  the  velocity  com¬ 
ponents. 

Fig.  21  shows  the  notation  used  in  9  ,  Z  plane.  From 
this  figure  it  can  be  seen  that  K  -  X  f  /\j  ,  where 

is  the  angular  distance  between  the  mean 
line  and  the  frozen  streamline  surface,  the  mean  line  lies 
half-way  between  two  adjacent  chord  surfaces.  Since  the 
blade  spacing  is  $  -  +  the  choice  of  K  and  ^ 

determines  yJ  ,  T"1**  and  6  • 

By  substituting  eqs.  (III. 1.9)  and  (III. 1.10)  into 
eqs.  (III. 1.1  to  III. 1.5)  after  eliminating  P  ,  ^  ,X 
by  using  eqs.  III.1.6)  and  (III.1.7).  Collecting  terms  of 
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equal  powers  of  S  and  equating  their  coefficients  to  zero,, 
one  obtains  five  simultaneous  algebraic  equations  for  each 
set  of  five  coefficients.  The  coefficients  of  5  up  to  the 
in  eqs.  (III. 1.1}  to  (III. 1.5)  were  calculated  and 
also  a  general  recurrence  formula  for  S71  with  ic 
The  solutions  of  these  equations  give  the  unknown  coeffi¬ 
cients  in  terms  of  ulq  ,  ,  /°e  ,  T©  ,  'LL,  , 

T,  and  either  v'[  or  u’t  .  The  three  arbitrary  func¬ 
tions  iLt  ,  ~T i  f  or  Zd',  may  be  used  to  ensure 

closure  of  the  blade  surfaces  at  the  trailing  edge. 

The  recurrence  equations  for  the  series  are  not  includ¬ 
ed  here  because  they  are  too  lengthy.  They  are  available  to 
interested  workers  at  the  Aeronautical  Engineering  Department, 
University  of  Illinois,  Urbana,  Illinois. 

III. 2.  TURBULENT  FLOWS 

The  general  Navier-Stokes  equations  were  used  with  the 
assumption  that  the  flow  is  composed  of  a  mean  flow  and  of 
turbulent  fluctuations  of  the  flow  variables.  This  gave  rise 
to  double  and  triple  correlations  between  these  variables. 

To  our  knowledge,  no  valid  theory  of  compressible  tur¬ 
bulent  flows  exists  which  would  allow  the  evaluation  of  these 
correlation  functions.  Therefore,  several  groups  in  this 
country  were  contacted  who  are  working  in  the  field  of  turbu¬ 
lence.  However,  it  was  pointed  out  by  these  workers  that  the 
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problem  cannot  be  solved  at  present  analytically*  Therefore, 


no  farther  work  was  carried  oat  on  this  phase* 
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IV.  CONCLUSIONS 


It  is  thought  that  the  report  presents  the  first  com¬ 
prehensive  treatment  of  turbomachine  analysis  including 
viscous  effects.  The  channel  theory  of  Lorenz  has  been 
developed  to  a  general  and,  it  is  hoped,  useful  stage. 

Both  compressible  and  incompressible  flows  are  treated 
which  may  be  laminar  or  turbulent  (except  that  the  general 
treatment  of  compressible  turbulent  flows  is  not  possible 
at  present). 

The  analyses  are  valid  for  axial  flow,  as  well  as 
radial  and  mixed  flow  turbomachines  (compressors  and 
turbines ) ,  although  the  numerical  work  included  treats 
axial  flow  compressor  stages  only.  This  is  so  because  the 
original  aim  of  the  work  was  the  development  of  a  theory 
for  axial  flow  turbomachines. 

It  is  realized  that  much  numerical  work  has  to  be 
done  to  render  the  analyses  useful  for  engineers.  Only 
after  such  work  has  been  carried  out  will  it  be  possible 
to  evaluate  the  usefulness  of  the  theory  presented. 
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TABLE  I.  DESIGN  PARAMETERS  FOR  NUMERICAL  EXAMPLE 


Length  of  stage 
Intake  radii 
Hub 

Intermediate 

Shroud 

Slope  of  shroud  (constant) 
Number  of  blades 
Parameter  in  eq#  (I.2.b.5) 
Pitch  angle 

Axial  velocity  component 

Intake  radial  velocity 
component  at  shroud 

Angular  velocity 

Eccentricity  of  velocity 
parabola 

Intake  pressure-density  ratios 
Hub 

Intermediate 

Shroud 

Density 


(4)  =  0,250  ft 

-  0.250  ft 

(fWt)  «  0.293  ft 

(tn)  =  0.333  ft 

(Ms)  *  0.5 
16 

(/9)  -  0.043 

(§)  *  tr/8  radians 

(xlTq)  »  100  ft/sec 

(uSL)  ■  50  ft/ sec 
Cco)  *  15,000  rpm 

(a)  -  0 

(PH.)«  862,770  ft2/sec2 
(PTOi)-  683,377  ft2/ sec2 
(PSi)  -  883,742  ft2/ sec2 
(/o)  s  0.002378  slugs/ft3 
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TABLE  II.  BLADE  TWIST  VARIATION  ALONG  AIIS 


z 

T 

0 

0 

.05 

.1274 

.10 

.1947 

.20 

.2712 

.30 

.3150 

.40 

.33^6 

.50 

.3622 

.60 

.3751 

.70 

.3636 

.BO 

.3690 

.90 

.3916 

1.00 

.3927 
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TABLES  III.  VARIATION  OF  RADIOS.  VELOCITIES.  AND 

- TOtsrogTOggf r  rato?  itorangH — * — 


Ilia.  Hub 


z/t 

r(ft ) 

u (ft/sec) 

v(ft/sec) 

V  (ft/sec) 

P(ft2/sec2) 

0 

.2500 

66.667 

392.700 

410.60 

082,770 

.05 

.2503 

65.741 

178.190 

214.64 

949,240 

.10 

.2664 

64.899 

113.250 

164.43 

964,040 

.20 

.2045 

63.420 

62.805 

134.04 

979,430 

.30 

.2902 

62.185 

41.096 

124.72 

991,800 

.40 

.3135 

61.122 

29.445 

120.04 

1,003,990 

.50 

.3207 

60.204 

20*444 

110.50 

1,016,200 

.60 

.3437 

59.404 

13.206 

117.07 

1,028,850 

.70 

.3505 

58.701 

8.623 

116.28 

1,041,730 

.00 

.3731 

50.080 

5.977 

115.80 

1,054,950 

.90 

.3075 

57.527 

2.775 

115.40 

1,060,560 

1.00 

.4010 

57.032 

0 

115.12 

1,082,530 

nib. 

Intermediate 

zA 

r(ft) 

u( ft/sec) 

v(ft/sec) 

V(ft/sec) 

P (ft2/ sec2) 

0 

.2933 

56.033 

460.650 

474.79 

803,380 

.05 

.3003 

56.533 

207.210 

236.92 

973,210 

.10 

.3074 

56.253 

130.660 

173.89 

991,440 

.20 

.3214 

55.747 

71.458 

134.96 

1,000,300 

.30 

.3353 

55.304 

46.209 

123.22 

1,021,060 

.40 

.3491 

54.912 

32.775 

118.70 

1,033,240 

.50 

.3627 

54.563 

22.557 

116.13 

1,045,550 

.60 

.3763 

54.252 

15.811 

114.86 

1,058,100 

.70 

.3899 

53.973 

10.631 

114.13 

1,070,970 

.00 

.4033 

53.721 

6.462 

113.70 

1,004,190 

.90 

.4167 

53.493 

2.984 

113.45 

1,097,780 

1.00 

.4301 

53.286 

0 

113.31 

1,111,730 

IIIc. 

Shroud 

zA 

r(ft) 

u(ft/sec J 

v(ft/sec) 

V (ft/sec) 

P (ft2/ sec2) 

0 

.3333 

50 

523.600 

535.40 

003,740 

.05 

.3396 

50 

234.280 

259.59 

998,560 

.10 

.3450 

50 

146.999 

184.69 

1,020,480 

.20 

.3503 

50 

79.672 

137.29 

1,038,970 

.30 

.3700 

50 

51.111 

122.93 

1,052,080 

.40 

.3033 

50 

35.994 

117.46 

1,064,360 

.50 

.3958 

50 

24.616 

114.48 

1,076,720 

.60 

.4083 

50 

17.155 

113.11 

1,089,200 

.70 

.4208 

50 

11.475 

112.39 

1,102,150 

.00 

.4333 

50 

6.943 

112.02 

1,115,360 

.90 

.4458 

50 

3.123 

111.85 

1,128,940 

1.00 

.4583 

50 

0 

111.80 

1,142,880 
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TABLE  IV.  VARIATION  OF  STATIC  AND  TOTAL 
PRESSURE-DENSITY  RATIOS  AT  EXIT 


re(ftj 

Pe(ft2/sec2) 

(ft2/ sec2) 
Ae 

Hub 

.4018 

1,082,530 

1,288,338 

Intermediate 

.4301 

1,111,730 

1,346,347 

Shroud 

.4583 

1,142,880 

1,408,287 

TABLE  V.  VARIATION  OF  STATIC  PRESSURE-DENSITY 
RATIO  AT  INTAKE 


rt(ft) 

Pi(ft2/sec2) 

Hub 

.2500 

882,770 

Intermediate 

.2933 

883,380 

Shroud 

.3333 

883,740 
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FIG.  I 

NOTATION  FOR  INFINITESIMAL  BLADE  SPACING 


FIG.  2 

RESULTANT  VELOCITY  DISTRIBUTION 
INFINITESIMAL  BLADE  SPACING 
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FIG.  4  4>,z  PROJECTION  OF  STREAMLINES 
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VELOCITY  IN  FT.  PER  SEC. 


FIG.  6  VELOCITY  COMPONENTS  VERSUS  f 
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FIG.  8  PRESSURE  RATIO  VERSUS  I 


-SHROUD 

-MEDIUM 

-HUB 


2  IN  FT. 


FIG.  II  TYPICAL  STREAMLINES  IN  THE  DEVELOPED 
HUB,  MEDIUM,  AND  SHROUD  SURFACES 
(BLADE  SHAPE) 
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SHROU 

MEDIUI 


VARIATION  OF  RESULTANT  VELOCITY  IN  THE 
MEDIUM,  AND  SHROUD  SURFACES 


-SHROUD 

-MEDIUM 


FIG.  13  AXIAL  VARIATION  OF  VELOCITY  COMPONENT! 
IN  THE  HUB,  MEDIUM,  AND  SHROUD 
SURFACES 
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FIG.  14  AXIAL  VARIATION  OF  STATIC  PRESSURE 

IN  THE  HUB,  MEDIUM,  AND  SHROUD  SURFACES 


SHROUD 

MEDIUM 


AXIAL  VARIATION  OF  DENSITY  IN  THE 
HUB,  MEDIUM,  AND  SHROUD  SURFACES 


- SHROUD 

- MEDIUM 

- HUB 


z  IN  FT. 


FIG.  16  TYPICAL  STREAMLINES  IN  THE  DEVELOPED 
HUB,  MEDIUM,  AND  SHROUD  SURFACES 
(BLADE  SHAPE) 

39 
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FIG.  17  AXIAL  VARIATION  OF  RESULTANT  VELOCITY  IN  THE 
HUB,  MEDIUM  AND  SHROUD  SURFACES 
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- MEDIUM 


HUB 


FIG.  18  AXIAL  VARIATION  OF  VELOCITY  COMPONENTS 
IN  THE  HUB,  MEDIUM,  AND  SHROUD 
SURFACES 
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FIG.  19  AXIAL  VARIATION  OF  STATIC  PRESSURE 

IN  THE  HUB,  MEDIUM,  AND  SHROUD  SURFACES 


VARIATION  OF  DENSITY  IN  THE 
MEDIUM,  AND  SHROUD  SURFACES 


:UPPER  CHORD  SURFACE 

LOWER  CHORD  SURFACE 

d>  :  FROZEN  STREAMLINE  SURFACE 
T0 

}jf  •  (r,z) 

FIG.  21  GEOMETRY  OF  BLADE  CHANNEL 
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